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ABSTRACT  (Continued) 

complex  beam  transport  system.  However,  radial  expansions  and  abrupt 
contraction  in  the  conducting  walls  plus  strong  axial  acceleration  in  the 
induction  gaps  of  the  accelerator  were  found  to  introduce  perturbations  to 
the  radial  beam  forces  which  drive  oscillations  in  the  beam  envelope.  The 
largest  radial  oscillations  were  excited  when  the  beam  was  injected  into  a 
strong  solenoidal  magnetic  field  to  compress  its  size  by  an  order  of  magnitude 
prior  to  transport  into  an  FEL  wiggler.  A  Hamiltonian  analysis  was  performed 
of  the  beam  cyclotron  oscillations  exciped  during  beam  transport.  It  was  shown 
that  the  beam  emittance  may  be  related  tc\  the  action  J  =  §  pdq,  which  is  pre¬ 
served  when  the  axial  gradients  in  the  radial  beam  forces  are  low  enough  that 
J  is  adiabatically  invariant.  However,  sudden  beam  compression  at  the  entrance 
to  the  wiggler  solenoid  was  predicted  to  excite  large  amplitude  cyclotron 
oscillations  on  the  beam  which  would  phase  mix  to  create  the  dominant  source  of 
emittance  growth  for  this  particular  transport  system.  Calculations  were  also 
performed  of  the  beam  quality  achievable  in  a  high  current,  magnetic  field- 
immersed  diode.  It  was  found  that  the  best  beam  quality  results  from  gently- 
shaped  anode  and  cathode  structures  which  create  electric  fields  whose  radial 
components  have  long  axial  gradient  lengths.  When  the  diode  axial  magnetic 
field  strength  is  high  enough  that  the  beam  electrons  are  tied  to  the  field 
lines  and  the  inverse  cyclotron  wave  number  is  below  the  axial  gradient  length 
of  the  radial  diode  forces,  extremely  high  quality  electron  beams  can  be  created. 
Beam  emittance  up  to  two  orders  of  magnitude  below  the  Lawson-Penner  emittance 
was  observed  in  some  cases. 
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I.  CALCULATIONS  OF  THE  FORMATION  AND  TRANSPORT  OF  HIGH 
QUALITY,  HIGH  CURRENT  RELATIVISTIC  ELECTRON  BEAMS 

Our  research  on  the  formation  and  transport  of  high 
quality,  high  current  relativistic  electron  beams  began  in 
August  1981  under  Contract  ONR  N00014-81-C-0704  and  continued 
for  four  years  until  July  31,  1985.  During  the  latter  portion 
of  this  contract,  work  was  also  initiated  on  the  subject  of 
the  phase  area  displacement  mode  of  free  electron  laser  opera¬ 
tion.  This  latter  work  has  been  continued  under  Contract  ONR 
N00014-86-C-0088  and  will  be  documented  in  the  reports  for  that 
contract.  Therefore,  this  final  technical  report  for  Contract 
ONR  N00014-81-C-0704  will  be  confined  to  the  subject  of  the 
formation  and  transport  of  high  quality,  high  current  electron 
beams. 

The  initial  impetus  for  research  on  this  topic  derived 
from  the  early  appreciation  of  the  difficult  tradeoff  which 
appeared  between  higher  current  beams  and  higher  quality  beams. 
For  a  variety  of  applications  of  relativistic  electron  beams, 
including  their  employment  as  drivers  of  free  electron  laser 
radiation,  performance  enhancement  requires  increasing  both  the 
beam  current  and  the  beam  quality.  Unfortunately,  it  was 
discovered  that  because  of  space  charge  effects,  higher  current 
beams  tend  to  possess  diminished  beam  quality  and,  conversely, 


higher  quality  beams  can  often  be  produced  only  at  the  expense  of 
some  reduction  in  beam  current. 


The  importance  of  the  electron  beam  driver  being  of  high 
quality,  with  low  emittance  and  low  energy  spread,  in  order  to 
achieve  high  gain  FEL  operation  and  reduce  competition  from  other 
radiation  mechanisms  driven  by  transverse  beam  thermal  energy  has 
long  been  well  recognized  (Refs.  1,2,3).  Although  higher 
electron  currents  tend  to  produce  higher  FEL  gain  due  to  more 
electrons  participating  in  the  FEL  interaction,  this  result  may 
not  be  achieved  if  the  beam  quality  is  reduced.  In  fact,  linear 
gain  can  plummet  if  the  electron  beam  quality  is  degraded  enough. 

The  tradeoff  which  must  be  confronted  between  beam 
current  and  beam  quality  is  illustrated  by  the  empirical 
Lawson-Penner  relation  (Refs.  4,5) 


£ 

—  =  Y  Sr  8  =  0.35 

TT 

where  £  is  the  normalized  beam  emittance.  Large  dividends 
n 

from  improved  FEL  performance  will  accrue  to  the  extent  that  this 
empirical  standard  can  be  surpassed  through  dedicated  and  clever 
beam  design.  It  was  one  of  the  primary  goals  of  the  research 
performed  under  this  contract  to  explore  the  physics  of  beam 
formation  and  transport  which  causes  the  beam  emittance  to  tend 
to  worsen  (i.e.,  increase)  with  increased  current.  It  was  hoped 
that  out  of  such  understanding  might  emerge  design  techniques  for 
delivering  relatively  lower  emittance  beams  at  higher  currents. 
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As  a  prototype  for  our  studies,  we  selected  the  Long 
Pulse  Induction  Linac  (LPIL),  which  was  built  at  the  National 
Bureau  of  Standards  (Ref.  6)  and  operated  at  the  Naval  Research 
Laboratory  (Refs.  3,7)  from  1980-1984.  The  device  is  illustrated 
in  Figure  1  and  had  nominal  operating  parameters  of  0.8  MeV, 

0.8  kA,  and  a  pulse  time  of  2  ysec.  The  electron  accelerator 
consisted  of  two  major  components,  an  injector  and  an  induction 
accelerator  module,  with  magnetic  focusing  coils  and  solenoids 
arrayed  axially  along  the  system  to  assist  in  beam  transport. 
Downstream  of  the  final  focusing  coil  shown  in  Figure  1,  the 
electron  beam  was  strongly  compressed  and  injected  into  a  wiggler 
solenoid  for  the  production  of  free  electron  laser  radiation. 

The  desired  objective  of  our  design  studies  was  to  learn 
how  to  best  form  a  high  current,  high  quality,  laminar  beam  in 
the  injector;  then  smoothly  transport  this  beam  several  meters 
downstream  through  the  induction  accelerator  module,  while 
preserving  as  much  of  the  beam  current  and  quality  as  possible; 
and  finally  compress  the  beam  strongly  and  inject  it  into  the 
wiggler  solenoid.  Unless  each  of  these  goals  are  well  achieved, 
the  fruits  of  success  in  the  early  stages  of  the  device  may  be 
lost  further  downstream  before  the  beam  delivers  its  payoff  in 
the  wiggler  interaction.  The  resultant  necessity  for  achieving  a 
comprehensive  design  for  beam  formation  and  transport,  which 
encompasses  the  entire  region  from  the  cathode  to  the  wiggler,  is 
common  to  a  wide  variety  of  similar  experimental  situations.  At 
the  same  time,  focusing  our  computational  and  theoretical  studies 


upon  this  specific  device,  from  which  experimental  data  was 
concurrently  available,  afforded  the  opportunity  to  achieve  those 
more  comprehensive  gains  in  understanding  which  flow  from  an 
active  interplay  between  theoretical  and  experimental  studies. 

The  results  of  our  studies  of  beam  formation  and 
transport  have  been  described  in  some  detail  in  Appendices  A 
through  H  hereto.  Appendices  A  through  D  comprise  four  letter 
progress  reports  which  describe  some  of  the  calculations  of  beam 
formation  in  the  diode,  as  well  as  many  of  the  detailed  calcula¬ 
tions  of  beam  transport  through  the  first  four  meters  of  the 
induction  linac.  Additional  calculations  of  beam  formation  and 
diode  transport  are  described  in  Appendix  E,  and  the  beam 
transport  through  the  accelerator  is  further  discussed  as  well. 

In  Appendix  F,  the  beam  formation  and  transport  is  summarized  and 
the  results  are  also  compared  with  the  predictions  of  a  one¬ 
dimensional  envelope  equation  for  R(z).  A  Hamiltonian  analysis 
is  also  presented  in  Appendix  F  of  the  excitation  of  cyclotron 
oscillations  upon  the  beam,  and  it  is  shown  how  one  may  strive  to 
design  the  transport  fields  to  minimize  these  cyclotron  waves. 
Appendix  G  contains  a  calculation  of  the  beam  emittance  which 
results  from  the  excitation  of  magnetic  cyclotron  oscillations 
during  beam  transport.  Finally,  Appendix  H  contains  a 
calculation  of  beam  quality  improvements  which  may  be  achieved 
through  the  use  of  high  magnetic  fields  and  gentle  wall  gradients 
in  field-immersed  diodes.  Below  we  will  mention  some  of  the 


highlights  of  the  research  described  in  these  Appendices  A 
through  H. 

Beam  formation  in  the  injector  was  studied  using  the 
REEFER  code  (Ref.  8)  and  the  STAGEN  code  (developed  at  Austin 
Research  Associates)  to  compute  steady  state  diode  flow  patterns 
in  r,z  geometry.  Early  examination  of  a  cantilevered  button 
field  emission  cathode  (shown  in  Figure  22  of  Appendix  A)  of 
14  cm  radius,  2.5  cm  thickness,  and  -400  kV  potential  inside  a 
grounded  cylindrical  housing  of  19  cm  radius  with  a  grounded 
anode  mesh  located  10  cm  downstream,  revealed  that  pre-breakdown 
electric  field  stresses  on  the  cathode  tip  would  exceed  170 
kV/cm.  This  was  found  to  lead  to  complete  breakdown  of  the 
cathode  surface,  a  current  density  profile  highly  peaked  on  edge, 
some  5  kA  of  emitted  current,  and  an  expanding  flow  pattern  in 
the  anode-cathode  gap  in  the  absence  of  a  strong  magnetic  guide 
field.  This  beam  was  judged  unacceptable  due  to  its  excessive 
size  (  >  14  cm  vs.  <  0.6  cm  desired  downstream  in  the  wiggler), 
its  excessive  current  (5  kA  vs.  —  1  kA  desired),  and  its 
expansive  flow  (see  Figure  23  of  Appendix  A)  in  the  diode — making 
laminar  compression  difficult  downstream.  Consequently,  the 
remainder  of  our  diode  studies  contemplated  the  use  of  either 
warm  or  cold  cathodes  with  annular  electrode  rings  in  the  anode- 
cathode  gap  to  provide  a  graded  potential  axially  through  the 
diode  and  to  reduce  the  pre-breakdown  stresses  on  the  cathode 
surface,  with  emission  out  to  —  8  cm  radius  controlled  by  the 
presence  of  fibers  on  the  cathode. 


Such  controlled  emission  cathodes  were  found  to  produce 


reasonable,  high  quality  beams  with  laminar  flow  patterns.  Diode 
parameters  selected  were 

anode-cathode  gap  spacing  =  d=  11.0  cm 
cathode  emission  radius  =  a  =  8.0  cm 
cathode  potential  =  -A<i>  =  -400  kv. 

The  gap  spacing  of  11  cm  is  sufficient  to  minimize  the  problem  of 
diode  closure  at  velocities  of  '■'*2.5  cm/usec  during  the  1-2  ysec 
pulse  time  of  the  beam.  The  emission  radius  of  8  cm  is  large 
enough  to  ensure  that  the  current  density  at  the  cathode  surface 
is  less  than  4  A/cm2,  as  required  for  reliable  operation  and  to 
avoid  the  Child-Langmuir  space  charge  limit  (see  calculation  in 
Appendix  E).  However,  the  requirement  that  the  beam  should 
eventually  be  compressed  to  a  radius  of  about  0.6  cm  in  the 
wiggler  solenoid  places  a  very  high  premium  on  techniques  for 
beam  convergence. 

It  is  easy  to  see  that  complete  reliance  on  magnetic 
compression  to  achieve  this  hundredfold  increase  in  the  beam 
current  density  is  hardly  feasible.  Hundreds  of  gauss  would  be 
required  to  tie  the  beam  well  to  field  lines  in  the  diode  region 
(see  magnetic  tying  formula  in  Appendix  E)  and  then  compression 
to  some  tens  of  kilogauss  would  be  needed  to  reduce  the  magnetic 
flux  surface  radius  to  0.6  cm  in  the  wiggler  solenoid. 
Consequently,  it  is  desirable  to  achieve  some  beam  convergence 
electrostatically  while  the  magnetic  field  is  weak,  leaving  a 

7 


lesser  amount  of  compression  to  be  achieved  magnetically. 
Electrostatic  beam  convergence  is  possible  in  a  graded  potential 
parallel  plate  geometry.  The  resistively  graded  electrodes  lead 
to  a  radial  boundary  condition  in  the  anode-cathode  gap  which  may 
be  modeled  as 


$  (r=b,  z)  -  <J> 


cathode 


A$ 


■  (i)E 
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By  varying  p,  one  may  vary  the  cathode  stresses,  the  radial 
current  density  profile,  and  the  diode  flow  pattern.  We 
evaluated  computationally  the  choices  p=l,  4/3,  1.52,  and  1.58, 
as  discussed  in  Appendices  C  and  E.  These  diode  potential 
gradings  are  illustrated  in  Figure  1  of  Appendix  C.  The  most 
desirable  results  were  obtained  for  a  sub-Child-Langmuir  grading 
of  p=1.52,  which  produced  relatively  weak  cathode  stresses, 
current  of  Ib  =  720  A,  and  a  convergent  flow  pattern  with  axially 
peaked  current  density  as  illustrated  in  Figures  2  and  3.  The 
convergent  beam  flow  occurs  because  the  graded  electrode 
structure  tends  to  short  out  the  radial  electric  fields  normally 
associated  with  the  beam  space  charge,  such  that  the  vz  x  B0 
pinch  force  is  left  to  compress  the  beam  radially. 

As  shown  in  Figure  1  (and  in  more  detail  in  Figure  1  of 
Appendix  A),  the  electron  beam,  once  formed,  must  be  transported 
through  a  complex  array  of  magnetic  field  coils  and  radially 
varying  conducting  boundaries,  and  through  two  accelerating  gaps 
in  the  induction  module  before  being  finally  compressed  and 


injected  into  the  wiggler  solenoid,  which  begins  380  cm 
downstream  of  the  cathode  surface.  In  order  to  study  this 
challenging  transport  problem,  we  employed  in  addition  to 
theoretical  analysis,  a  one-dimensional  R(z)  beam  envelope  code 
and  the  two-dimensional  (r,z),  steady  state,  relativistic 
particle  simulation  code  STAGEN. 

The  STAGEN  code  retains  all  three  components  of  velocity 
and  of  the  electric  and  magnetic  fields,  but  is  restricted  to  the 
two  independent  variables  r,z,  appropriate  to  situations  which  • 
are  symmetric  in  the  azimuthal  coordinate  6.  in  order  to  reduce 
the  computing  requirements  to  a  manageable  level,  the  axial 
transport  is  divided  into  four  "stages."  The  first  stage  extends 
to  z  *  23  cm,  encompassing  the  11  cm  length  of  the  anode-cathode 
region  plus  12  cm  beyond  the  anode  mesh.  The  flow  is  initiated 
from  an  equipotential  near  but  not  on  the  cathode  surface — 

typically  ^ _<{>cathode  -  kV*  sta9e  2  extends  from  z  -  23  cm 
to  z  »  111  cm,  encompassing  the  initial  focusing  coil  and  the 
conductor  wall  flaring  out  to  a  radius  of  17.5  cm.  Stage  3 
extends  from  z  *  111  cm  to  z  *  258  cm,  and  encompasses  three 
more  focusing  coils,  the  first  induction  accelerating  gap  (of 
200  kV),  and  the  first  induction  solenoid.  Stage  4  extends  from 
z  *  258  cm  to  z  *  406  cm,  and  encompasses  the  second  induction 
solenoid,  the  second  induction  accelerating  gap  (also  of  200  kV), 
another  focusing  coil,  and  the  first  26  cm  of  the  wiggler 


solenoid. 


Under  this  contract,  STAGEN  was  exercised  repeatedly  to 
examine  the  beam  transport  through  the  first  four  meters  of  the 
long  pulse  induction  linac  for  a  variety  of  diode  potential 
gradings  and  for  several  magnetic  field  configurations.  Because 
the  constraints  of  operating  with  the  existing  coil  configuration 
and  achieving  an  order  of  magnitude  beam  radius  compression  into 
the  wiggler  solenoid  make  the  transport  so  difficult,  emphasis  on 
the  STAGEN  runs  was  first  placed  on  simulating  typical  experi¬ 
mental  parameters.  In  this  way,  support  for  the  experiment  was 
provided,  feedback  was  received,  and  an  appreciation  for  the 
crucial  elements  of  the  transport  physics  was  obtained.  An 
example  of  these  transport  simulations  is  provided  in  Figures  4 
through  7.  Figure  4  displays  the  magnetic  field  profile  in  the 
LPIL  for  a  typical  set  of  coil  currents,  including  the  magnetic 
ramp  into  the  2  kG  wiggler  solenoid.  A  separate  magnetic  field 
code  was  written  and  exercised  in  order  to  generate  magnetic 
field  profiles  such  as  that  shown  in  Figure  4,  from  the  known 
coil  currents,  dimensions,  and  locations  (shown  in  Figure  1  of 
Appendix  A ) . 

Figure  5  shows  the  equipotentials  in  the  diode  region  and 
in  the  12  cm  beyond  the  anode  mesh.  It  may  be  seen  that  the 
graded  electrode  structure  does  an  excellent  job  of  shorting  out 
the  radial  beam  space  charge  electric  field  in  the  diode  region, 
allowing  the  gently  convergent,  laminar  beam  flow  seen  in 
Figure  2.  A  theoretical  calculation  was  made  of  the  rms  scatter 
which  would  occur  in  the  transverse  beam  velocity  as  the  beam 
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Figure  6.  Transport  and  Compression  of  LPIL  Beam 
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passed  through  the  anode  mesh  at  z  *  11  cm,  due  to  the  charge 
induced  on  the  wire  mesh.  It  was  estimated  that 


1  rms 


1  b 

2/i"  d 


(3) 


where  b  is  the  distance  between  wires  in  the  anode  mesh,  d  is 

the  anode-cathode  separation,  and  (y  -  1)  me2  is  the  diode 

injection  energy.  For  diode  parameters  of  400  keV  injection 

energy,  with  one  rail  tungsten  anode  mesh  wires  on  a  grid  of  98% 

transparency,  one  finds  that  b/d  =  0.025  and  ySB  ___  *  0.005, 

x  r  ins 

corresponding  to  emittance  growth  of  about  ~25ir  mrad-cm  at  the 
anode  mesh.  However,  this  anode  mesh  emittance  growth  is  small 
in  comparison  with  the  experimentally  estimated  emittance  of 
~200tt  mrad-cm  created  at  the  surface  of  the  hairbrush  cathodes. 

To  mock  up  the  effect  of  this  cathode-generated  emittance,  a 
scatter  corresponding  to  200tt  mrad-cm  emittance  was  therefore 
introduced  at  the  anode  mesh  for  these  runs.  Figure  6  displays 
the  beam  transport  through  Stages  1-4  of  the  LPIL  device,  and 
the  beam  compression  and  injection  into  the  wiggler  solenoid. 
(Full-scale  versions  of  these  figures  appear  as  Figures  5  and  7 
through  9  of  Appendix  F. )  It  may  be  seen  that  135A  of  current  is 
lost  to  the  wall  just  beyond  the  first  induction  gap.  For  this 
particular  case,  the  second  induction  gap  was  inactive.  The 
strong  beam  compression  of  the  remaining  583A  of  current  is  seen 
in  Stage  4  (Figure  9  of  Appendix  F).  However,  this  strong 
compression  excites  large  amplitude  zero  frequency  cyclotron 
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oscillations  which  are  displayed  in  the  $t~z  phase  space  plot 

of  Figure  7.  This  strong  sausaging  will  subsequently  phase  mix 

to  create  a  downstream  emittance  of  ~940tt  mrad-cm  for  this 

transport  example.  In  Figure  8,  the  variation  of  y  with  z 

z 

is  displayed  for  the  electrons  in  Stage  4.  At  the  throat  of  the 

wiggler,  the  median  value  of  y  is  found  to  be  «1.6j  there  is 

z 

a  considerable  spread  since  the  outer  electrons  have  been  excited 
to  much  higher  transverse  velocities  than  have  the  inner 
electrons.  These  details  on  the  phase  space  distribution  of  the 
electrons  turned  out  to  be  quite  useful  (Ref.  3)  in  interpreting 
that  certain  of  the  experimental  radiation  measurements  indeed 
signify  an  FEL  mode,  as  expected. 

Examination  of  these  STAGEN  calculations  of  beam 
transport  over  several  meters  distance  with  multiple  magnetic 
coils,  multiple  accelerating  gaps, and  variable  wall  geometry 
(appropriate  to  the  LPIL  accelerator  at  NRL)  revealed  several 
sources  of  beam  quality  degradation.  Noteworthy  is  the 
significant  wall  flaring  at  z  =  42-60  cm,  the  abrupt  reduction 
of  wall  radius  near  z  *  160  cm,  the  strong  axial  electric  fields 
applied  at  the  two  induction  gaps,  and  the  very  strong  magnetic 
compression  at  the  entrance  to  the  wiggler  solenoid.  Each  of 
these  transport  features  can  perturb  the  radial  force  balance  of 
the  electron  beam,  so  as  to  excite  radial  beam  oscillations  which 
may  lead  to  emittance  growth  (and  in  some  case  to  current  loss  as 
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In  order  to  theoretically  analyze  emittance  growth  in  a 
complex  transport  scenario  such  as  this,  we  found  it  useful  to 
employ  a  Hamiltonian  treatment  (Ref.  9)  of  the  axially-dependent , 
radial  confining  fields.  This  Hamiltonian  may  be  made  consistent 
with  a  one-dimensional  envelope  equation  for  R(z),  similar  to 
the  equation  developed  by  Lee  and  Cooper  (Ref.  10).  This 
analysis  is  described  in  Appendix  F  and  provides  guidance  for 
improving  the  beam  transport.  The  assumptions,  equations,  and 
symbols  used  in  this  analysis  are  summarized  in  Tables  1  and  2  of 
Appendix  F.  A  normalized  form  of  the  envelope  equation  appears 

as 
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and  where  the  beam  energy  |y(z)  -  lj  me2  is  obtained  from 
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where 


(Y  -  1)  me 2  =  diode  injection  energy 


b(z)  =  conducting  wall  radius 
wall 


e$.._-,i(z)  *  wall  potential. 


The  primary  source  of  beam  quality  degradation  during 
beam  transport  is  the  excitation  of  beam  cyclotron  oscillations 
as  a  result  of  (frequently)  abrupt  axial  variations  in  the  radial 
beam  focusing  forces.  The  beam  equilibrium  radius  is  that  for 


which  the  radial  "force"  terms  of  Equation  (4)  are  in  balance, 


when  <R"  =  — j  (R  =  0.  The  physical  significance  of  these  force 
dz 


terms  is  discussed  in  Table  1  (extracted  from  Appendix  F).  The 
Y'2  term  is  seen  to  represent  primarily  the  focusing  kicks  that 
the  beam  receives  in  the  diode  and  in  the  induction  gaps.  During 


the  early  portion  of  the  transport,  while  Bz  £  200  G,  the  beam 


radius  is  determined  by  competition  between  the  magnetic  V^BZ 


focusing  force  and  the  defocusing  space  charge  force  [terms  (1) 
and  (3)  of  Table  1] .  However,  as  the  magnetic  field  is  increased 
toward  2  kG,  the  defocusing  space  charge  force  is  outweighed  by 
the  centrifugal  and  pressure  forces  [term  (4)  of  Table  1].  Then 
the  beam  radius  is  determined  by  competition  between  this  term 
(4)  and  the  magnetic  pinch  term  (1)  of  Table  1,  leading  to  a 
smaller  equilibrium  radius  at  very  large  magnetic  fields. 
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Table  1 


RELATIVE  SIGNIFICANCE  OF  TERMS  IN  ENVELOPE  EQUATION 


(1) 


Y28  V 


—  Magnetic  v,pBz  focusing  forces  which  are 
important  everywhere  except  in  the  diode 
region  where  Bz  is  very  weak. 


(2) 
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—  Radial  focusing  forces  which  are  significant 
(for  beam  current  <<  Alfven  current)  only  in 
the  diode  and  in  the  induction  gaps.  These 
forces  represent  the  replacement  of  Er  by 
Ez  fields,  and  tend  to  locally  unbalance  the 
radial  force  equilibrium. 


(3) 
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Defocusing  space  charge  forces  which  are  in 
competition  with  magnetic  focusing  forces 
during  the  first  three  meters  of  large  radius 
beam  transport. 


c  +  P 

(4)  i>  —  Centrifugal  and  pressure  forces  which  are  in 

competition  with  magnetic  focusing  forces 
during  the  strong  magnetic  compression  of  the 
beam  to  a  very  small  radius.  Both  c  and 
are  constant  downstream  of  the  anode,  and 
e diode  =  200-300  mrad-cm  dominates  the 
P<*  =  50  mrad-cm  due  to  the  weak  leakage 
magnetic  field  at  the  cathode. 
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In  order  to  transport  a  fairly  cool  laminar  beam,  one 
would  desire  the  focusing  forces  to  vary  rather  slowly  in  z,  so 
that  the  beam  never  received  large  radial  kicks  corresponding  to 
large  departure  from  radial  force  balance.  However,  many  beam 
transport  systems,  including  the  LPIL  system,  do  not  afford  the 
luxury  of  gentle  axial  variations.  For  example,  in  the  injector 
diode  and  particularly  in  the  induction  gaps,  large  V  *  dY/dz 
appears  and  tends  to  locally  imbalance  the  radial  force 
equilibrium.  Likewise,  when  the  conductor  walls  abruptly  flare 
or  contract  (as  they  may  be  seen  to  do  in  Figure  6),  this  forces 
the  beam  to  propagate  across  equipotentials,  which  introduces  y ’ 
via  the  logarithmic  space  charge  term  of  Equation  (5).  When  the 
beam  transport  is  magnetically  focused,  axial  variations  in  the 
magnetic  field  strength  will  directly  unbalance  the  force  equi¬ 
librium.  When  force  equilibrium  is  unbalanced,  the  beam  will 
subsequently  oscillate  about  the  equilibrium  radius.  For 
magnetically  focused  transport,  these  oscillations  will  occur  at 
the  cyclotron  frequency.  Due  to  radial  variations  in  the 
particle  energy,  these  cyclotron  oscillations  will  subsequently 
phase  mix  to  produce  an  effective  increase  in  the  beam 
temperature,  or  transverse  emittance. 

The  principal  emittance  growth  observed  during  the 
transport  of  the  beam  through  the  LPIL  system  occurred  due  to  the 
strong  excitation  of  cyclotron  oscillations  during  magnetic 
compression  at  the  mouth  of  the  wiggler  solenoid.  The  magnetic 
compression  is  shown  in  Figure  4,  the  tight  beam  cyclotron 
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oscillations  may  be  seen  in  Figure  6,  near  z  =  400  cm,  and  the 
corresponding  large  perturbations  in  the  radial  electron  velocity 
may  be  seen  in  Figure  7.  The  beam  was  also  observed  to  undergo 
significant  radial  oscillations  and  occasional  tight  pinching  at 
earlier  stages  of  the  transport.  For  example,  Figure  6  displays 
a  fairly  tight  beam  pinch  near  z  *  140  cm  under  the  third 
focusing  coil.  For  some  transport  cases  computed  under  this 
contract,  the  pinching  was  so  severe  that  a  space  charge  virtual 
cathode  was  established  and  some  electrons  were  reflected 
backward  from  the  pinch  to  limit  the  current  transmitted.  Such 
"limiting  current"  phenomena  can  occur  when  the  logarithmic  space 
charge  term  of  Equation  (5)  becomes  so  large  that  electrons  are 
slowed  and  Equation  (5)  cannot  be  obeyed  unless  the  beam  current 
is  reduced.  Radial  beam  oscillation  and  pinching  can  also  result 
in  the  loss  of  some  beam  current  to  the  walls,  as  illustrated  in 
Figure  6  near  z  *  200  cm. 

We  have  also  performed  a  Hamiltonian  analysis  of  the 
cyclotron  oscillations  excited  upon  the  electron  beam,  which  is 
summarized  in  Table  4  of  Appendix  F.  It  is  shown  that  the  beam 
envelope  oscillations  may  be  viewed  as  anharmonic  oscillations  in 
the  pseudopotential  P(R),  illustrated  in  Figure  12  of 
Appendix  F.  One  may  analytically  evaluate  the  effect  of  sudden 
or  adiabatic  variations  in  the  force  coefficients  of  the 
Hamiltonian,  so  as  to  assess  the  magnitude  of  transverse 
oscillations  produced  on  the  beam  during  transport.  The  beam 
emittance  may  be  related  to  the  action  J  =  ^p  dq,  and  is  found 


to  be  preserved  when  changes  occur  slowly  so  that  J  is  also 
adiabatically  invariant. 

On  the  other  hand,  the  magnetic  compression  illustrated 
near  z  =  400  cm  in  Figure  4  occurs  suddenly  rather  than 
adiabatically,  and  was  observed  to  produce  large  cyclotron 
oscillations  (see  Figure  7).  Theoretically,  it  is  indeed  found 
that  large  cyclotron  oscillations  of  the  magnitude  observed 
should  be  created  if  the  magnetic  field  is  suddenly  increased 
greatly  at  a  point  where  the  beam  radius  is  relatively  large. 
Appendix  G  contains  a  calculation  of  the  beam  emittance  which 
results  from  the  phase  mixing  of  cyclotron  oscillations  excited 
by  the  propagation  of  a  beam  through  a  sudden  magnetic 
discontinuity. 

Some  experimental  attempts  were  made  to  comply  with  the 
theoretical  recommendations,  by  introducing  an  additional 
transport  coil  to  make  the  magnetic  compression  more  adiabatic — 
and  indeed  some  improvement  in  transport  was  seen.  However,  the 
continued  necessity  for  substantial  magnetic  compression  in  order 
to  inject  into  the  small  wiggler  solenoid,  a  large  portion  of  the 
beam  created  in  a  large  area,  low  magnetic  field  diode,  implied  a 
certain  irreducible  amount  of  cyclotron  oscillation  emittance 
growth  at  the  mouth  of  the  wiggler.  (The  constraint  of  working 
with  the  existing  accelerator  configuration  presented  a  limita¬ 
tion  to  the  improvements  which  could  be  practically  realized  by 
adding  or  modifying  the  transport  field  coils.)  These  induced 
transverse  cyclotron  oscillations  were  in  turn  found 


experimentally  to  produce  a  high  level  of  cyclotron  radiation, 
while  the  emittance  growth  reduced  the  level  of  FEL  radiation 
which  could  be  obtained.  Consequently,  a  decision  was  eventually 
taken  by  NRL  experimentalists  to  eliminate  the  strong  magnetic 
compression  into  the  wiggler  (which  hence  became  free  of  the 
guide  field),  thereby  aperturmg  awa /  most  of  the  transported 
beam.  This  step  essentially  eliminated  the  unwanted  cyclotron 
radiation,  although  the  FEL  radiation  amplitude  was  consequently 
low  because  of  the  low  driving  current. 

The  remainder  of  our  research  under  this  contract  on  the 
production  and  transport  of  high  quality,  high  current  electron 
beams  was  concentrated  on  a  calculation  of  the  beam  quality 
achievable  in  high  current,  magnetic  field-immersed  diodes.  This 
calculation  is  presented  in  Appendix  H.  We  found  that  the  best 
beam  quality  results  from  gently-shaped  anode  and  cathode 
structures,  which  create  electric  fields  whose  radial  components 
have  long  axial  gradient  lengths.  When  the  diode  axial  magnetic 
field  strength  is  high  enough  that  the  beam  electrons  are  tied  to 
the  field  lines  and  the  inverse  cyclotron  wave  number  8zc/^  is 
below  the  axial  gradient  length  of  the  radial  diode  forces, 
extremely  high  quality  electron  beams  can  be  created. 

Improvements  in  the  beam  emittance  of  one  to  two  orders  of 
magnitude  over  the  Lawson-Penner  emittance  of  Equation  (1)  were 
observed  in  some  cases.  In  one  calculation  of  diode  emittance 
from  a  smoothly  tapered  anode  structure,  it  was  found  that 
doubling  the  magnetic  guide  field  caused  the  transverse  electron 


velocity  to  decrease  by  a  factor  of  nine — which  illustrates  the 
potential  of  magnetic  field-immersed  diodes  for  producing  high 
quality  electron  beams. 


We  are  hopeful  that  the  research  completed  under  this 
contract  has  confirmed  the  importance  of  enlightened  beam 
transport  design  to  the  goal  of  producing  high  quality,  high 
current  beams  and,  further,  that  it  has  indicated  certain  avenues 
toward  the  achievement  of  this  goal  through  the  use  of  gentle 
axial  variation  in  the  transport  structure,  and  strong  focusing 
fields. 
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Or.  C.  W.  Roberson 

Plasma  Physics  Division,  Code  4704 

Naval  Research  Laboratory 

4555  Overlook  Avenue 

Washington,  D.C.  20375 

Dear  Chuck: 

Enclosed  is  some  computer  output  to  illustrate  some  of 
our  preliminary  results  on  the  beam  transport  problem  through  the 
first  four  meters  of  the  NBS  induction  accelerator,  and  also  some 
preliminary  diode  calculations  with  the  tentative  cold  cathode 
design  which  you  supplied.  I  have  numbered  the  figures  and  will 
provide  some  brief  commentary  keyed  to  those  numbers. 

Figure  1  shows  the  layout  of  the  machine.  The  cathode 
has  a  radius  of  about  8  cm  and  is  located  at  z  *  0.  The  wiggler 
solenoid  begins  at  z  *  380  cm.  We  have  moved  the  break  point 
between  Stages  1  and  2  to  z  *  23  cm,  as  shown.  Stage  4  extends 
from  z  ■  258  cm  to  z  *  406  cm,  which  is  26  cm  into  the  wiggler 
solenoid. 


Figure  2  shows  the  magnetic  field  profile  which  we  used 
for  our  initial  transport  studies.  It  does  not  include  the 
strong  2  kG  magnetic  field  in  the  final  solenoid.  Figure  3  shows 
how  the  magnetic  field  profile  is  changed  in  the  presence  of  the 
2  kG  solenoidal  field.  However,  we  have  only  just  today  sub¬ 
mitted  our  first  runs  to  attempt  to  propagate  the  beam  through 
the  strong  magnetic  compression  in  the  wiggler  solenoid,  and  we 
don't  have  the  results  back  yet. 

Figure  4  shows  the  Stage  1  propagation  from  near  the 
cathode  through  the  anode  screen  and  into  the  first  drift 
section.  The  electrons  are  initiated  at  20  keV  energy,  1.2  cm 
downstream  of  the  true  cathode  surface,  in  accordance  with  the 
Child-Langmuir  law.  The  vertical  tic  marks  represent  the 
resistively  graded  electrodes  between  the  cathode  and  anode.  The 
electrode  potentials  were  in  accordance  with  the  z1*  Child- 
Langmuir  law,  to  yield  vertical  equipotentials  for  space-charge 
limited  flow.  A  current  of  850  amperes  was  emitted,  and  the  A-K 
potential  jump  was  400  kV.  Scatter  was  not  introduced  in  the 
initial  electron  momenta  to  mock  up  the  anode  screen  scatter  or 
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cold-cathode  emission  effects.  A  uniform  current  density  profile 
was  assumed.  Hence  this  run  might  represent  a  hot  cathode  run, 
had  you  been  able  to  extract  850  A  from  your  hot  cathode. 

Figure  5  shows  the  equipotential  contours.  It  is  clear 
that  the  graded  electrode  structure  does  a  good  job  in  eliminat¬ 
ing  Er  within  the  A-K  gap.  This  is  why  the  beam  flow  in 
Figure  4  is  initially  convergent,  since  the  vz  B^  pinch  force 
is  dominant  in  the  A-K  gap.  Near  z  ■  20  cm,  the  usual  Er  beam 
fields  have  been  established. 

Figure  6  shows  the  magnetic  flux  contours.  The  leakage 
field  in  the  A-K  gap  region  is  around  5-10  G. 

Figures  7  and  8  show  8r  and  8q  for  the  Stage  1  flow. 
Note  that  the  sign  should  be  reversed  on  all  of  the  vertical  axis 
quantities  that  are  plotted. 

Figure  9  shows  the  Stage  2  flow  pattern.  The  Stage  2 
z-coordinates  are  relative  to  the  start  of  Stage  2  at  z  *  23  cm 
(i.e.,  z  -  23  cm  should  be  the  label).  The  zero  frequency  oscil¬ 
lations  which  are  set  up  as  the  beam  propagates  into  the  18  cm 
radius  tube  and  encounters  the  first  magnetic  field  peak  near 
z  =  72  cm  (z  -  23  cm  =  49  cm)  are  clearly  evident. 

Figure  10  shows  the  Stage  2  flux  lines,  while  Figures  11 
and  12  show  the  8r  and  8q  flow  field.  the  beam  seems  to  be 
transporting  fairly  well  here,  without  too  much  phase  mixing.  At 
the  end  of  Stage  2,  the  counter-current  in  coil  2  has  driven  the 
B z  field  near  zero,  so  that  8q  is  back  near  zero. 

Figure  13  shows  the  Stage  3  flow  pattern.  During  this 
stage,  the  electron  beam  propagates  through  two  Bz  peaks  and 
across  the  first  200  kV  induction  gap  at  z  -  182  cm  (or 
z  -  111  cm  =  71  cm) .  The  outer  beam  electrons  come  very  near 
the  wall  at  z  =  154  cm  and  z  *  258  cm.  Figure  14  shows  the 
equipotential  contours,  including  the  effect  of  the  induction 
gap.  Figure  15  shows  the  Stage  3  flux  contours. 

Figures  16  and  17  show  3r  and  Bg  for  Stage  3.  The 
beginnings  of  phase  mixing  are  here  apparent. 

Figure  18  shows  the  Stage  4  flow  pattern,  encompassing 
the  second  induction  gap,  and  the  320  G  magnetic  field  peak  under 
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coil  7.  The  beam  expands  near  the  end  of  Stage  4  due  to  the 
absence  of  the  strong  wiggler  solenoidal  magnetic  field. 

Figure  19  shows  the  magnetic  flux  surfaces,  revealing  the  fall- 
off  in  Bz  near  z  *  400  cm,  just  as  in  Figure  2. 

Figures  20  and  21  show  8r  and  8@  in  Stage  4,  and 
considerable  phase  mixing  seems  to  be  apparent  here.  The  beam 
"thermal  velocity"  component  may  be  a  few  %c  here.  It  remains  to 
be  seen  how  well  this  beam  can  be  compressed  to  ~2  kG  by  using 
the  Figure  3  magnetic  profile,  and  what  its  emittance  then  is. 

Subsequent  transport  runs  will  include  introducing 
initial  beam  scatter  to  gauge  the  transport  sensitivity  to  this, 
studying  transport  through  the  strong  magnetic  compression  into 
the  wiggler,  and  trying  to  numerically  repeat  the  Bz  profile 
optimization  experiments  to  some  degree. 

The  remaining  figures  concern  our  initial  code  diode 
studies.  We  began  by  considering  the  tentative  design  which  you 
had  given  us,  shown  in  Figure  22.  The  12.7  cm  dimension  was 
changed  to  14  cm,  upon  the  assumption  that  the  figure  was  drawn 
to  scale.  The  guide  magnetic  field  was  set  to  zero,  and  an  anode 
screen  was  placed  10  cm  downstream  of  the  cathode  face. 

Figure  23  shows  the  preliminary  flow  field  obtained  from 
the  REEFER  code.  Because  of  the  absence  of  electrode  grating  or 
of  pierce  electrode  structure,  the  flow  is  somewhat  expansive. 

We  have  arbitrarily  cutoff  emission  at  r  =  11  cm,  in  agreement 
with  the  outermost  fibers  shown  in  Figure  22.  However,  one 
problem  which  is  apparent  is  that  field  stresses  on  the  cathode 
tip  are  measured  to  be  about  170  kV/cm,  and  it  would  seem  likely 
that  the  cathode  would  field  emit  all  the  way  out  to  the  tip 
under  these  stresses.  This  is  in  rough  agreement  with  analytic 
estimates  of  ~200  kV/cm  at  the  tip.  The  stresses  are  this  high 
because  of  the  thinness  ( Az  is  small)  of  the  cathode  tip,  and  its 
nearness  (4.7  cm)  to  the  inside  surface  of  the  outer  wall.  The 
cathode  potential  is  assumed  to  be  -400  kv  w.r.t.  ground  at  the 
outer  wall  and  at  the  anode  screen. 

Because  of  the  high  electric  field  stresses  near  the 
outer  portion  of  the  cathode,  the  current  pulled  off  in  an  edge- 
peaked  profile,  as  shown  in  Figure  24.  Across  most  of  the 
cathode  face,  the  current  density  is  about  5  to  6  A/cm2,  in  crude 
agreement  with  the  Child-Langmuir  expectation  for  parallel  plate 
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geometry.  All  total,  3.3  kA  are  pulled  out  of  the  cathode.  It 
is  tempting  to  consider  aperturing  the  beam  down  to  an  8  cm 
radius,  thus  transmitting  the  central  1.2  kA  or  so,  which  has  a 
much  more  uniform  current  profile. 

Figure  25  shows  <J>  (  r=0  ,  z  )  /V0  plotted  against  z/d  ,  where 
V0  =  400  kV  and  d  =  10  cm.  Apparently,  <j>  <x  z1,28,  in  close 
agreement  with  the  z ^  of  Child-Langmuir . 

Subsequent  runs  will  include  introducing  a  magnetic  guide 
field  into  the  diode  region,  exploring  foilless  situations,  and 
considering  other  cathode  designs. 


3 

a 


We  realize  that  we  are  still  in  the  preliminary  stages 
of  this  study  of  beam  transport  and  diode  design,  but  we  are 
optimistic  that  the  pace  of  obtaining  physics  results  may  quicken 
somewhat  now.  I  hope  these  results  will  nevertheless  be 
interesting  and  helpful  to  you,  and  we  will  keep  you  informed  of 
our  further  progress. 


Sincerely, 


/f7  - 


Bob  Thompson 


P.S.  The  ARA  participants  in  this  study  are  Barry  Moore, 
John  Uglum,  Lee  Sloan,  and  myself. 
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Figure  2.  Magnetic  Field  Profile  Thru  Accelerator 
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Figure  19.  Magnetic  Flux  Contours  in  Stage  4 
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Letter  from  Dr.  James  R.  Thompson  to  Dr.  C. 
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April  27,  1982 


Dr.  C.  W.  Roberson 

Plasma  Physics  Division,  Code  4704 

Naval  Research  Laboratory 

4555  Overlook  Avenue 

Washington,  D.C.  20375 

Dear  Chuck: 

Enclosed  is  some  more  computer  output  to  illustrate  some 
of  our  recent  calculations  on  beam  transport  through  the  first 
four  meters  of  the  NBS  induction  accelerator,  including  injec¬ 
tion  into  the  high  field  wiggler  solenoid.  I  have  numbered  the 
figures  and  the  following  commentary  is  keyed  to  those  numbers. 

Figure  1  shows  the  magnetic  field  profile  which  we  used 
for  the  transport  calculations  included  herein.  The  wiggler 
magnetic  field  was  assumed  to  be  2  kG,  and  at  the  end  of  our 
Stage  4,  at  z  »  406  cm  downstream  of  the  cathode,  Bz  has  risen 
to  about  1.96  kG.  This  point  is  26  cm  into  the  throat  of  the 
wiggler  solenoid. 

Figures  2  through  5  show  the  flow  patterns  for  a  cool, 
laminar  beam  as  it  propagates  through  our  Stages  1,  2,  3,  and  4 
in  the  presence  of  the  magnetic  field  of  Figure  1.  Note  that  the 
axial  scale  varies  from  stage  to  stage,  and  the  computer  numbers 
are  referenced  to  the  beginning  of  each  stage.  In  Figure  2,  the 
v2  Ba  pinching  in  the  diode  is  apparent,  due  to  Er  being 
mostly  shorted  out  there.  At  the  end  of  Stage  3,  in  the  center 
of  the  two  induction  modules,  some  electrons  graze  the  wall  and 
subsequently  collapse  to  small  radii  in  Stage  4.  Near  the  end  of 
Stage  4,  in  the  region  of  strong  magnetic  compression,  some 
electrons  are  spun  up  to  high  8-velocities,  and  about  72  A  of 
current  is  lost  on  the  wall.  However,  about  780  A,  or  92%  of  the 
current  was  successfully  transmitted  into  the  solenoid. 

Figures  6  and  7  show  8r  and  69  in  the  Stage  4  flow. 
Strong  zero  frequency  cyclotron  oscillations  are  apparent  in  the 
region  of  the  high  field  ramp.  At  this  point,  most  of  the 
transverse  beam  motion  is  still  coherent,  or  sloshing,  but  some 
thermalization  has  clearly  begun.  If  one  assumes  the  worst,  that 
all  of  this  cyclotron  motion  will  eventually  thermalize,  then  the 
resulting  beam  emittance  would  be  as  follows: 
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This  is  considerably  higher  than  the  Stage  1  emittance  for  the 
laminar  flow  shown  in  Figure  2.  There  we  estimate 
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The  reason  for  the  much  higher  emittance  in  Stage  4  is  of  course 
the  strong  zero  frequency  oscillations  which  are  excited  as  the 
beam  is  compressed. 

Figures  8  and  9  show  the  equipotentials  and  the  magnetic 
flux  surfaces  in  Stage  4. 

Figure  10  shows  the  Stage  1  diode  flow,  with  a  random 
angular  scatter  introduced  at  the  location  of  the  anode  mesh. 

The  RMS  angular  scatter  was  taken  as  2.83°  in  the  transverse 
plane.  Figures  11  and  12  show  Sr  and  Bq  for  this  scattered 
Stage  1  flow.  The  emittance  just  downstream  of  the  anode  mesh  is 
as  follows: 
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Figure  13  shows  the  equipotentials  in  the  diode  region. 

Figure  14  shows  the  Stage  2  flow  pattern  of  the  scattered 
beam,  and  Figures  15  and  16  show  the  corresponding  3r  and  B0 . 

Figure  17  shows  the  Stage  3  flow  pattern,  and  Figures  18 
and  19  show  the  corresponding  Sr  and  Bq  .  Stage  3  encompasses 
coil  2,  coil  3,  coil  4,  the  first  induction  gap,  and  the 
solenoidal  coil  5.  In  this  stage,  a  significant  portion  of  the 
scattered  beam  is  lost  to  the  walls.  About  135  A,  or  16%  of  the 
beam  current  is  lost  at  z  =  154  cm  at  the  entrance  to  the 
induction  module.  Another  349  A,  or  41%  of  the  beam  current  is 
lost  between  20  cm  and  70  cm  downstream  of  the  first  induction 
gap,  under  the  solenoidal  coil  5.  This  leaves  374  A,  or  44%  of 
the  beam  current  injected  into  stage  4.  Figure  20  shows  the 
Stage  3  equipotentials,  including  the  strong  gradients  in  the 
induction  gap. 

Figure  21  shows  the  Stage  4  flow  pattern  for  the 
scattered  beam,  and  Figures  22  and  23  show  the  corresponding  Br 
and  Bq.  Stage  4  encompasses  the  second  induction  gap,  and  the 
magnetic  field  ramp  into  the  wiggler  solenoid.  No  further 
current  is  lost  in  Stage  4.  The  Br  and  Bq  diagnostics  again 
reveal  strong  zero  frequency  cyclotron  oscillations.  These  have 
begun  to  thermalize  somewhat,  but  are  still  largely  coherent  at 
the  end  of  Stage  4.  If  we  again  assume  that  this  cyclotron 
motion  will  eventually  phase  mix,  the  resulting  emittance  would 
be  as  follows: 


Dr.  C.  W.  Roberson 
Naval  Research  Laboratory 
April  27,  1982 
Page  4 


5  8  RMS 
a 
R 

Y 

£n 


0.20  I 

1.10  cm  j 

0.78  cm 
2.57  (800  keV) 

400  mrad-cm  @  I  = 


eyeball  estimates 


374  A 


The  projected  emittance  for  the  scattered  beam  is  thus 
the  same  as  for  the  unscattered  beam;  however,  since  the 
transmitted  current  is  much  lower  for  the  scattered  beam,  the 
scattered  beam  may  be  said  to  be  hotter  than  is  a  corresponding 
amount  of  current  in  the  unscattered  beam.  Likewise,  for  the 
374  A  transmitted,  the  emittance  has  grown  significantly  between 
Stages  1  and  4  due  to  the  cyclotror.  motion  excited  in  Stage  4. 

Figure  24  shows  the  magnetic  flux  contours  for  Stage  4. 
The  inner  contour  drawn  corresponds  to  Bzr  -  0.72  kG-cm2 . 

Hence  it  corresponds  to  the  f .x  wit -in  r  *  0.6  cm  at 
B  z  =  2  kG,  which  are  the  des.:«*  ‘  . -.al -state  beam  conditions. 
Unfortunately,  Figure  21  revea..-  'tat  the  beam  is  well  outside  of 
this  desired  flux  contour,  3esp.-e  arrying  less  current  than 
desired  into  Stage  4.  Producina  the  desired  beam  current  density 
in  the  wiggler  solenoid  with  b.~  =  J  nG  therefore  remains  a 
difficult  challenge  for  our  futjre  work. 


I  am  sorry  for  the  delay  of  several  days  in  getting 
these  results  into  the  mail.  Also,  we  have  recently  noticed 
peculiarity  in  the  current  diagnostic  which  appears  to  have 
relation  to  the  current  lost  in  stage  3.  until  we  get  this 
sorted  out,  you  should  probably  treat  the  figure  of  374  A  of 
transmitted  current  with  a  grain  of  salt. 


a 

some 


Sincerely, 


Bob  Thompson 
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Figure  3.  Laminar  Beam  Flow  a  Stage  2 
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Figure  8.  Equipotentials  for  Compression  of 
Laminar  Beam  a  Stage  4 


Flux  Contours  for  Stage  4  Compression 
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Figure  15.  Scattered  Flow  3  Stage  2 
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Figure  17.  Scattered  Flow  Thru  Stage  3,  With  Current  Loss 
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Figure  18.  Scattered  Flow  Thru  Stage  3 
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Figure  21.  Transmission  of  Reside 
Current  Thru  Stage  4  0 


Compression  of  Scattered  Beam  Flow  a  Stage  4 
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APPENDIX  C 


Letter  from  Dr.  James  R.  Thompson  to  Dr.  C. 
dated  Hay  13,  1982. 


includes  30  figures 
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Dr.  C.  W.  Roberson 
1908  Tollbridge  Court 
Alexandria,  Virginia  22308 

Dear  Chuck: 

Enclosed  is  the  latest  batch  of  copies  of  our  computer 
output,  to  illustrate  our  recent  progress  in  calculating  diode 
flow  and  beam  transport  through  the  first  four  meters  of  the  MBS 
induction  accelerator.  1  have  numbered  these  figures  and  the 
following  commentary  is  keyed  to  those  numbers. 

We  have  done  a  series  of  diode  calculations,  to  determine 
the  effect  of  various  alternative  potential  gradings  of  the  ring 
electrodes  in  the  anode-cathode  gap  on  the  beam  formation  and  on 
the  subsequent  transport.  Figure  1  illustrates  the  four 
different  gradings  which  we  have  considered.  The  "design"  and 
the  "modified  design"  cases  correspond  to  the  two  sets  of  numbers 
on  the  engineering  drawing  of  the  electrode  geometry.  Notice 
that  the  linear  grading  produces  stronger  electric  fields  near 
the  cathode  than  does  Child-Langmuir  grading,  while  the  design 
grading  produces  weaker  electric  fields.  As  expected,  the 
extracted  current  is  ordered  in  the  same  way  as  the  electric 
field  stresses. 

Figure  2  shows  the  diode  flow  for  linear  grading;  1.24  kA 
current  is  extracted.  Figure  3  shows  the  cathode  current  density 
profile  for  this  case.  Because  of  the  high  electric  field 
stresses,  the  current  density  is  peaked  at  the  edge  of  the  beam. 
The  diode  flow  is  also  slightly  divergent  near  the  beam  edge. 

Figure  4  reviews  the  r-z  transport  geometry,  and  shows 
the  locations  of  our  computer  stage  break-points.  Figure  5 
reviews  the  magnetic  field  profile  which  is  based  upon  the  data 
that  you  gave  us  in  December,  1981,  and  which  we  have  been  using 
until  very  recently. 

Figures  6  through  9  illustrate  the  transport  of  the  beam 
formed  in  a  linearly  graded  diode  (see  Figures  2  and  3),  and 
propagated  along  the  system  shown  in  Figures  4  and  5.  No  anode 
mesh  scattering  was  introduced.  Out  of  1,260  amperes  initially 
extracted  from  the  cathode,  720  amperes  are  lost  in  Stage  3  just 
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beyond  the  first  induction  gap.  The  remaining  540  amperes  (43% 
of  the  current)  are  then  well  compressed  in  Stage  4. 

Figures  10  and  11  show  the  transverse  velocities  in  the 
Stage  4  compression  to  3Z  =  2  kG.  Assuming  as  usual  that  the 
cyclotron  oscillations  will  eventually  phase  mix,  the  resulting 
emittance  would  be  as  follows: 


RMS 


a 


0.07 
0.52  cm 


R  =  3l/J~2  -  0.37  cm 


y  «  2.57  (800  keV) 

n 

—  =66  mrad-cm  @  I  =  540  A 

IT 

This  residual  beam  is  of  high  quality,  and  the  Figure  9  flow 
pattern  may  be  seen  to  be  compressed  within  the  desired  inner 
flux  contour  of  Figure  12. 


Figure  13  shows  the  diode  flow  pattern  for  a 
Child-Langmuir  potential  grading.  The  beam  is  slightly 
convergent,  and  the  current  density  profile  is  only  slightly 
peaked  on  edge,  as  shown  in  Figure  14.  Figures  15  through  22 
review  previous  transport  calculations  for  Child-Langmuir 
potential  grading.  However,  for  these  previous  transport 
studies,  the  cathode  current  density  profile  was  assumed  uniform 
in  radius:  the  slight  peak  shown  in  Figure  14  was  omitted. 
Likewise,  the  current  was  only  about  850  A  rather  than  900  A  as 
in  Figures  13  and  14. 

Figures  15  through  18  show  the  transport  without  anode 
mesh  scattering.  figures  19  through  22  show  the  same  transport 
with  2.83°  of  RMS  scattering  at  the  anode  mesh.  This  is  a  fairly 
large  amount  of  mesh  scattering.  Figure  21  shows  that  135  A  were 
lost  at  the  entrance  to  the  first  induction  solenoid.  Although 
this  case  does  not  correspond  closely  with  the  experimental 
cases,  this  location  is  where  you  have  apparently  been  losing  a 
lot  of  current  also  in  your  experimental  runs.  An  additional 
349  A  were  lost  downstream  of  the  first  induction  gap,  where  you 
have  also  been  losing  a  smaller  amount  of  current  experimentally. 
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Figure  23  shows  the  diode  flow  for  the  design  potential 
grading,  which  is  a  sub-Child-Langmuir  grading.  For  this  case, 
we  allowed  the  cathode  to  emit  all  the  way  out  to  the  tip,  to  see 
how  much  additional  current  might  be  produced.  We  got  820  A, 
with  a  current  density  profile  peaked  on  axis  as  shown  in 
Figure  24. 

Figure  25  shows  the  diode  flow  for  the  set  of  parameters 
that  you  recently  gave  us.  We  used  the  modified  potential 
grading,  and  the  most  up-to-date  geometrical  dimensions.  The 
cathode  emission  extended  to  r  =  a  =  8.0  cm.  The  leakage  mag¬ 
netic  fields  from  coil  #1  were  5.7  gauss  at  the  cathode  and 
8.6  gauss  at  the  anode  mesh.  The  flow  is  fairly  strongly 
convergent,  and  720  A  were  extracted,  which  is  some  10%  lower 
than  you  observe.  Figure  26  shows  that  the  cathode  current 
density  is  peaked  on  axis  and  declines  to  about  85%  of  the  peak 
value  at  the  edge  of  the  beam. 

We  have  just  started  a  transport  calculation  for  the 
beam.  We  are  using  the  magnetic  field  profile  shown  in 
Figure  27,  based  upon  the  coil  settings  which  we  obtained  from 
you  during  our  February,  1982  visit.  We  are  again  assuming  2  k'j 
in  the  wiggler  solenoid. 

Figure  28  shows  the  Stage  1  diode  flow  for  this  beam.  We 
have  introduced  anode  mesh  scattering  of  0rras  ~  0.28°,  which  is 
much  less  than  for  our  previous  scattering  runs,  but  which  we 
believe  may  be  more  realistic.  The  modified  design  for  the  ',ub- 
Child-Langmuir  potential  grading  has  been  used.  Figures  29  and 
30  show  the  perpendicular  velocities  for  the  Stage  1  flow.  At 
this  point  the  beam  flow  is  still  extremely  laminar. 

Our  computer  is  presently  chewing  on  the  transpjrt 
through  Stages  2,  3,  and  4.  I  hope  to  have  some  of  these  results 
to  send  to  you  tomorrow,  or  at  least  to  relate  over  the 
telephone.  I'll  also  try  to  get  busy  on  our  renewal  proposal. 
Hope  that  you  have  a  good  meeting. 
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Figure  3.  Current  Density  Profile  a  Cathode 
for  Linear  Potential  Grading 


Laminar  Beam  Formation  a  Stage  1,  with  Linear 
Potential  Grading  [o  to  400  kvL  Edge-Peaked  Jz(r) 
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Figure  8.  Transport  Thru  Stage  3,  with  Current  Loss  to  Walls 
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Figure  12,  Flux  Contours  for  Stage  4  Compression 
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Figure  14.  Current  Density  Profile  a  Cathode 

for  Child-Langmuir  Potential  Grading 
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Figure  21.  Scattered  Flow  Thru  Stage  3,  with  Current  Loss 
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Figure  27.  Magnetic  Field  Profile  from  Feb.  1982  Data 


7-W-8i ? 

//://:£?  122 
Cor- 


STAGE  1 

fSITION 

.  OOO  1 6 . 000  20 . 000  24 . 000 


STAGE  1 

Z-POSITION 

.000  4.000  8.000  12.00016.00020.00024.000 


f 


M 


I 

I 

I 

§ 

1 

I 

fi 

1 

I 

I 

B 

1 

i 

i 

I 


BETA  0 


Figure  30.  Stage  1  with  Mesh  Scattering  eDMc  =  0.28 
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APPENDIX  D 


Letter  from  Dr.  Barry  N.  Moore  to  Dr.  C.  W.  Roberson, 
dated  June  29,  1982. 


§ 


—  includes  25  figures 
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Dr.  C.  W.  Roberson 

Plasma  Physics  Division,  Code  4704 

Naval  Research  Laboratory 

4555  Overlook  Avenue 

Washington,  D.C.  20375 

Dear  Chuck: 

Enclosed  is  computer  output  from  two  more  series  of 
calculations  of  beam  transport  through  the  first  four 
meters  of  the  long  pulse  induction  linac.  I  have  numbered 
these  figures  and  the  following  commentary  is  keyed  to  those 
numbers. 

The  first  series  of  runs,  in  Figures  1-14,  simulates 
transport  of  a  cold  cathode  beam,  with  1.4°  of  rms  angular 
scattering  introduced  at  the  anode  mesh  to  mock  up  the 
emittance  generated  at  the  hairbrush  cathode  face.  Figure  1 
shows  the  magnetic  field  profile  which  was  used  for  these 
transport  studies.  Figure  2  shows  the  computed  beam  flow 
in  the  diode,  using  the  experimental  sub-Child-Langmuir 
potential  grading.  We  find  that  720  amperes  of  current  is 
extracted,  peaked  slightly  on  axis  as  shown  in  Figure  3. 

The  current  density  near  the  beam  edge  is  about  85%  of  the 
central  current  density.  Figure  4  shows  the  convergent 
beam  flow  through  Stage  1,  including  the  anode  mesh 
scattering.  Figure  5  shows  the  transport  through  Stage  2, 
and  Figure  6  shows  the  beam  focusing  and  injection  into  the 
first  induction  module  in  Stage  3.  It  may  be  seen  that 
about  135  amperes  or  19%  of  the  beam  current  is  lost  to  the 
walls  just  downstream  of  the  first  induction  gap.  In 
Figure  7,  the  beam  transport  through  Stage  4  is  shown,  for 
a  case  where  the  second  induction  module  has  been  turned 
off.  The  strong  compression  of  the  residual  current  by  the 
2  kG  wiggler  magnetic  guide  field  is  evident  near  the  end 
of  Stage  4.  Figures  8  and  9  show  the  transverse  velocity 
modulations  of  the  beam  particles  in  Stage  4,  while 
Figure  10  shows  the  Yz  factors  for  the  beam  electrons  in 
Stage  4.  The  strong  cyclotron  oscillations  of  the  beam 
near  the  end  of  Stage  4,  at  the  entrance  of  the  wiggler 
solenoid,  may  be  seen  to  extract  energy  from  the  directed 
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z-motion  and  to  convert  it  to  r,  0  motion.  It  appears 
that  the  rms  beam  radius,  averaged  over  the  cyclotron  wave, 
is  about  R  =  1.2  cm.  The  rms  value  of  8x,  averaged  over 
the  cyclotron  wave,  appears  to  be  8X  rms  ~  0.39.  Since 
Y  »  2,  the  normalized  emittance  is  about 

~T  ~  YRBx  rms  ~  940  mrad-cm 

after  the  cyclotron  wave  phase  mixes.  The  residual  beam 
current  for  this  case  is  583  amperes.  It  may  be  noticed 
that  there  is  a  considerable  spread  in  the  Yz  values  at 
the  entrance  to  the  wiggler,  since  the  outer  electrons  have 
considerably  more  transverse  cyclotron  motion  than  the 
inner  electrons  do.  The  median  value  of  Yz  appears  to  be 
about  1.61  at  the  entrance  to  the  wiggler  for  this  case. 

In  Figures  11-14,  the  Stage  4  transport  was 
recomputed  with  the  second  induction  gap  turned  on.  The 
higher  beam  energy  created  may  be  seen  to  lengthen  the 
cyclotron  wavelength,  as  expected.  In  this  case,  it 
appears  that  R  =  1.3  cm,  8X  rms  =  °-4°r  Y  =  2.4, 

(en/Tr)  =  1260  mrad-cm  and  Yz_me(jian  =  1.76. 

In  Figures  15-25,  we  have  repeated  a  transport 
calculation  for  a  high  current,  warm  cathode  beam. 

However,  we  applied  the  December  1981  scaling  factors  to 
the  coil  currents  rather  than  to  the  fields  themselves, 
yielding  the  magnetic  field  profile  shown  in  Figure  15. 

The  magnetic  field  goes  slightly  negative  between  105  cm 
and  125  cm  from  the  cathode,  reaching  perhaps  -32G  under 
the  second  focusing  coil.  This  implies  two  weak  magnetic 
cusps  located  at  105  cm  and  125  cm. 

Figure  16  shows  the  hot  cathode  beam  formation  in 
Stage  1;  it  was  assumed  that  858A  were  emitted.  The 
experimental  sub-Child-Langmuir  potential  grading  was  used, 
with  uniform  current  density  profile.  Figure  17  shows  the 
beam  transport  in  Stage  2,  through  the  first  magnetic  cusp 
at  z  -  105  cm.  As  the  beam  approaches  the  magnetic  cusp, 
its  radius  increases  somewhat,  but  its  azimuthal  velocity 
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returns  to  zero  and  goes  slightly  negative  as  shown  in 
Figure  18.  The  azimuthal  velocity  is  low  near  the  cusp, 
since  the  beam  was  born  in  a  weak  field  region.  Because  of 
the  large  beam  radius  acquired  in  Stage  2,  the  beam  cannot 
be  sufficiently  compressed  in  Stage  3  to  avoid  the  loss  of 
significant  current  at  the  entrance  to  the  first  induction 
module.  In  Figure  19  it  may  be  seen  that  600A  or  70%  of 
the  beam  current  is  lost,  and  only  258A  or  30%  of  the 
current  is  injected  through  the  induction  module.  Figure 
20  shows  the  electron  Yz  factors  in  Stage  3;  the  loss  of 
the  outer,  transversely  excited  electrons,  and  the 
acceleration  of  the  central  core  of  the  beam  is  apparent. 
Figure  21  shows  the  transport  of  this  residual  beam  current 
through  Stage  4,  and  the  beam  compression  at  the  entrance 
to  the  wiggler  solenoid.  Figures  22  and  23  show  the 
transverse  beam  velocity  components  in  Stage  4,  while 
Figures  24  and  25  show  8Z  and  Yz  in  Stage  4.  In 
this  case,  the  cyclotron  wave  amplitude  is  much  less  than 
before,  but  there  is  still  an  appreciable  spread  in  Yz 
due  to  the  cyclotron  motion.  From  Figures  21-25,  we 
estimate 


R  =  0.6  cm 


6 


x  rms 


0.21 


—  =  310  mrad-cm 


^z-median 

Of  course,  the  transmitted  current  of  258A  is  much 
lower  for  this  case  than  for  the  cold  cathode  case  discussed 
earlier . 
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We  are  proceeding  with  envelope  equation  studies 
and  with  studies  of  emittance  generation  by  the  cold, 
hairbrush  cathodes.  Of  course,  it  should  be  recognized 
that  the  envelope  equations  will  be  of  limited  usefulness, 
since  they  cannot  represent  the  effects  of  radial 
gradients,  such  as  in  Y(r).  We  will  keep  you  informed  of 
our  progress. 


Sincerely, 

Barry  NO«oore 


BNMrkft 


Enclosures:  As  noted  above 


Figure  1.  Magnetic  Field  Profile  from  Feb.  1982  Data 
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Figure  15,  Magnetic  Field  Profile  from  Dec.  1981  Data 
with  Scaled  Coil  Currents 
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Figure  16,  Hot  Cathode  Beam  Formation  in  Stage  1,  with  Sub-Chi ld-Langmuir  Potential 
Grading  (Modified  Design),  and  Uniform  Current  Density  Profile 


Figure  18.  Laminar  Beam  Rotational  Velocity  in  Stage  2,  as 
Beam  Passes  Through  a  Vfeak  Magnetic  Cusp 
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igure  24.  Acceleration  and  Compression 
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Figure  25.  Acceleration  and  Compression  of  Residual  Beam  Flow  in  Stage  4 


APPENDIX 


Review  of  Beam  Formation  and  Transport  Through  an  Induction  Linac 


Section  II  of  Renewal  Proposal  No.  I-ARA-82-U-56  (ARA-454), 
"Computational  and  Theoretical  Studies  of  the  Formation  and 
Transport  of  High  Quality,  High  Current  Electron  Beams  for  Use 
as  Free  Electron  Laser  Drivers,"  May  1982. 
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REVIEW  OF  BEAM  FORMATION  AND  TRANSPORT 
THROUGH  AN  INDUCTION  LINAC 

Results  of  our  calculations  of  beam  formation  and 
transport  through  the  system  illustrated  in  Figure  1  have 
previously  been  discussed  during  one  site  trip  to  the  Naval 
Research  Laboratory  and  in  three  letter  progress  reports,1 
and  will  be  briefly  reviewed  herein. 

Beam  formation  in  the  injector  has  been  studied  by 
2 

using  the  REEFER  code  to  compute  steady  state  diode  flow 
patterns.  One  design  with  a  cantilevered  button  field 
emission  cathode  of  14  cm  radius,  2.5  cm  thickness,  and 
-400  kV  potential  inside  a  grounded  cylindrical  housing  of 
19  cm  radius  with  a  grounded  anode  mesh  located  10  cm  down¬ 
stream  was  evaluated.  However,  the  pre-breakdown  electric 
field  stresses  were  found  to  be  above  170  kV/cm  on  the 
cathode  tip,  which  would  lead  to  complete  breakdown  of  the 
cathode  surface,  a  current  density  profile  highly  peaked  on 
edge,  some  5  kA  of  emitted  current,  and  an  expanding  flow 
pattern  in  the  anode-cathode  gap  in  the  absence  of  a  strong 
magnetic  guide  field.  This  beam  was  judged  unacceptable, 
inasmuch  as  its  size  {>14  cm)  was  so  much  greater  than  the 
desired  downstream  beam  radius  in  the  wiggler  (<0.6  cm),  th 
flow  was  expansive  in  the  diode — making  laminar  compression 
difficult  downstream,  and  the  current  was  well  above  the 


nominal  design  value  of  %  1  kA.  Consequently,  the  remainder 
of  our  diode  studies  have  contemplated  the  use  of  either 
warm  or  cold  cathodes,  with  emission  out  to  ^ 8  cm  radius 
controlled  by  the  presence  of  fibers,  and  with  annular 
electrode  rings  in  the  anode-cathode  gap  to  provide  a  graded 
potential  axially  through  the  diode,  and  to  reduce  the 
pre-breakdown  stresses  on  the  cathode  surface. 

The  following  general  theoretical  considerations 
reveal  the  principal  constraints  which  govern  the  diode  and 
beam  transport  design.  First,  the  phenomenon  of  diode 
closure,  at  velocities  of  ^  2  cm/ysec  during  the  1-2  ysec 
pulse  time  of  the  beam,  leads  to  the  necessity  for  an 
anode-cathode  gap  spacing  d  of  at  least  10  cm  to  avoid 
excessive  changes  in  the  beam  impedance. 

Second,  it  is  necessary  that  the  diode  operate  at  a 
current  density  below  the  Child-Langmuir  limit.  For  a 
"parallel-plate"  diode,  which  the  resistively  graded 
electrode  rings  simulate  to  some  extent,  one  has 

I  me  3/e 

Tta2  ~  87rd2 

where  me  3/e  =  17  kA.  For  a  400  kV  injector, 
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and  one  finds  that  the  beam  current  density  J  must  not 
exceed  about  5  A/cm2  for  d  >  10  cm.  Therefore,  in  order  to 
produce  1,000  amperes  of  current,  the  emitting  radius  of  the 
cathode  must  exceed  about  8  cm. 

Third,  the  requirement  that  the  electron  beam  should 
eventually  be  compressed  to  a  radius  of  about  0.6  cm  in  the 
wiggler  solenoid  places  a  very  high  premium  on  techniques 
for  beam  convergence.  It  is  easy  to  see  that  complete 
reliance  on  magnetic  compression  to  achieve  this  hundredfold 
increase  in  the  beam  current  density  is  hardly  feasible.  In 
order  for  the  beam  to  be  magnetically  compressed,  it  must 
first  be  well  tied  to  the  field  lines,  which  roughly  occurs 
when 
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where  gj  ,  ft  are  the  relativistic  plasma  and  cyclotron 

r 

frequencies.  This  tying  requirement  would  require  hundreds 
of  gauss  in  order  to  be  obeyed  in  the  diode  region.  Then 
if  the  beam  were  to  be  magnetically  compressed  by  an  order 
of  magnitude  in  radius  while  tied  to  the  flux  surface 
Bza2  =  constant,  the  magnetic  field  would  have  to  be 
increased  to  some  tens  of  kilogauss  in  the  wiggler  solenoid. 
Judging  this  to  be  technologically  infeasible  or  too 


expensive,  we  conclude  that  the  beam  convergence  must  be 
achieved  electrostatically  in  part,  while  the  magnetic  field 
is  weak,  leaving  a  lesser  amount  of  compression  to  be 
achieved  magnetically. 

Electrostatic  beam  convergence  is  possible  either 
through  the  use  of  Pierce  electrode  geometry  or  in  a  graded 
potential  parallel  plate  geometry.  We  have  concentrated  on 
this  latter  possibility,  which  is  built  into  the  existing 
long  pulse  induction  linac.  Diode  parameters  for  this 
injector,  which  we  have  adopted  for  our  computational 
studies  as  well,  are 

anode-cathode  gap  spacing  =  d  =  11.0  cm 
cathode  emission  radius  =  a  =  8.0  cm 
cathode  potential  =  -  =  -400  kV. 

The  resistively  graded  electrodes  lead  to  a  radial  boundary 
condition  in  the  anode-cathode  gap  which  may  be  modeled  as 


$  (r=b,  z)  -  4> 


cathode 


We  have  considered  computationally  three  alternative  cases: 


(i)  p  =  1,  linear  grading 

(ii)  p  =  4/3,  Child-Langmuir  grading 

(iii)  ps:  1.52,  sub-Child-Langmuir  experi¬ 
mental  grading. 

The  computed  diode  flow  patterns  and  the  current  density 
profiles  at  the  cathode  surface  are  illustrated  for  these 
three  cases  in  Figures  2  through  7.  A  grounded  anode  mesh 
was  placed  at  z  =  d  =  11  cm,  to  conform  to  the  experimental 
situation.  The  axial  magnetic  field  was  varied  in  the  range 
0-10  G,  commensurate  with  the  magnitude  of  the  initial 
experimental  magnetic  fields,  and  the  diode  flow  was  found 
to  be  insensitive  to  Bz  in  this  range.  Just  as  the 
electric  field  stresses  in  the  vicinity  of  the  cathode  are 
strongest  for  p  =  1  and  weakest  for  p  =  1.52,  the  current 
extracted  was  found  to  be  highest  for  p  =  1  (I  =  1.24  kA) 
and  lowest  for  p  =  1.52  (I  3  0.72  kA) .  It  is  also  found 
that  for  the  linear  grading  (p  =  1),  the  current  density  was 
strongly  peaked  at  the  edge  of  the  beam  and  the  outer  flow 
was  somewhat  divergent.  For  the  intermediate  Child-Langmuir 
grading  (p  =  4/3) ,  the  current  density  is  only  slightly 
peaked  at  the  edge  of  the  beam  and  the  flow  is  slightly 
convergent.  (The  computed  current  of  I  =  0.90  kA  is  also 
in  close  agreement  with  the  theoretical  Child-Langmuir 
prediction  of  0.89  kA.)  Finally,  for  the  sub-Child-Langmuir 
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grading  (p  -  1.52),  the  current  density  is  slightly  peaked 
on  axis  and  the  flow  pattern  is  strongly  convergent,  as 
illustrated  in  Figure  6.  The  convergent  beam  flow  results 
from  the  graded  electrode  structure  tending  to  short  out 
the  radial  electric  fields  normally  associated  with  the 
electron  beam  space  charge  such  that  the  vz  x  Bq  pinch 
force  is  left  to  compress  the  beam  radially. 

As  shown  in  Figure  1,  the  electron  beam,  once 
formed,  must  be  transported  through  a  complex  array  of  mag¬ 
netic  field  coils  and  radially  varying  conducting  boundaries, 
and  through  two  accelerating  gaps  in  the  induction  module 
before  being  finally  compressed  and  injected  into  the 
wiggler  solenoid,  which  begins  380  cm  downstream  of  the 
cathode  surface. 

In  order  to  study  this  challenging  transport 
problem,  Austin  Research  Associates  has  further  extended  and 
developed  a  steady  state,  fully  relativistic  particle  simu¬ 
lation  code  named  STAGEN.  This  code  retains  all  three 
components  of  velocity  and  of  the  electric  and  magnetic 
fields,  but  is  restricted  to  the  two  spatial  coordinates 
r,z  as  independent  variables,  appropriate  to  situations 
which  are  symmetric  in  the  azimuthal  coordinate  9.  In 
order  to  reduce  the  computing  requirements  to  a  manageable 
level,  the  axial  transport  is  divided  into  four  "stages." 

The  first  stage  extends  to  z  =  23  cm,  encompassing  the 


11  cm  length  of  the  anode-cathode  region  plus  12  cm  beyond 
the  anode  mesh.  The  flow  is  initiated  from  an  equi- 
potential  near  but  not  on  the  cathode  surface — typically 
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0-<J>  .  <20  kV.  Stage  2  extends  from  z  =  23  cm  to 

z  =  111  cm,  encompassing  the  initial  focusing  coil  and  the 
conductor  wall  flaring  out  to  a  radius  of  17.5  cm.  Stage  3 
extends  from  z  =  111  cm  to  z  -  258  cm,  and  encompasses 
three  more  focusing  coils,  the  first  induction  accelerating 
gap  (of  200kV) ,  and  the  first  induction  solenoid.  Stage  4 
extends  from  z  =  258  cm  to  z  -  406  cm,  and  encompasses 
the  second  induction  solenoid,  the  second  induction 
accelerating  gap  (also  of  200  kV) ,  another  focusing  coil, 
and  the  first  26  cm  of  the  wiggler  solenoid. 

Under  this  contract,  STAGEN  has  been  exercised 
repeatedly  to  examine  the  beam  transport  through  the  first 
four  meters  of  the  long  pulse  induction  linac  for  all  three 
of  the  diode  potential  gradings  discussed  previously,  and 
for  several  magnetic  field  configurations.  We  have  also 
developed  the  capability  of  introducing  scattering  into  the 
beam  flow  as  it  passes  through  the  anode  mesh.  To  illustrate 
these  transport  calculations.  Figure  8  displays  one  magnetic 
field  configuration  based  upon  some  December  1981  coil 
settings,  and  including  the  magnetic  ramp  into  the  2  kG 
wiggler  solenoid.  For  this  magnetic  field  configuration  and 
for  the  Child-Langmuir  potential  grading  of  Figure  4,  the 
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beam  transport  is  shown  in  Figures  9  through  12  for  Stages  1 
through  4.  Figure  9  shows  the  slightly  convergent  diode 
flow,  with  rms  mesh  scattering  of  2.83°  introduced.  A 
current  of  858  A  was  followed  in  this  series  of  runs,  which 
is  slightly  below  the  900  A  given  by  the  diode  calculations. 
This  current  is  transported  into  Stage  3,  but  there  135  A 
(or  15.7%)  of  the  current  is  lost  at  the  entrance  to  the 
induction  module,  while  an  additional  349  A  (or  40.7%)  of 
the  current  is  lost  downstream  of  the  first  induction  gap, 
as  shown  in  Figure  11.  The  remaining  374  A  (or  43.6%)  of 
the  current  is  transported  through  Stage  4  and  injected  into 
the  wiggler  solenoid  without  further  loss,  as  shown  in 
Figure  12.  The  r  and  e  velocity  flow  patterns  of  the 
compressed  Stage  4  flow  are  shown  in  Figures  13  and  14.  It 
is  apparent  that  strong  zero  frequency  cyclotron  oscilla¬ 
tions  have  been  excited  in  the  magnetic  compression  to  2  kG. 
These  cyclotron  oscillations  have  begun  to  phase  mix  and 
will  continue  to  thermalize  further  downstream  in  the 
wiggler.  From  these  flow  diagnostics,  the  eventual  beam 
emittance  may  be  estimated  as  follows: 

RMS  transverse  velocity/c  =  66^  -  0.20 
RMS  beam  radius  =  R  =0.78  cm 
Y  =  2.57  (800  keV) 

£ 
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This  illustrated  transport  calculation  is  neither  the  most 
careful  attempt  to  duplicate  experimental  conditions  nor  an 
optimized  transport  design.  However,  it  does  illustrate  the 
capabilities  of  the  STAGEN  code.  Experimental  measurements 
have  been  reported  which  indicate  beam  current  loss  at  the 
same  two  locations  indicated  in  Figure  11,  but  we  have  not 
yet  duplicated  the  experimental  situation  well  enough  to  be 
able  to  claim  good  quantitative  agreement.  Such  transport 
calculations  are  currently  underway,  using  the  sub- 
Child-Langmuir  potential  grading  which  is  experimentally 
in  place  and  which  yields  more  convergent  flow  in  the 
injector  diode. 
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Review  of  Previous  Results  on  Beam  Transport  Studies 
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REVIEW  OF  PREVIOUS  RESULTS 


ON  BEAM  TRANSPORT  STUDIES 

Results  of  our  calculations  of  beam  formation  and 
transport  through  the  system  illustrated  in  Figure  1  have 
previously  been  discussed  during  several  site  trips  to  NRL,  in 
several  letter  progress  reports,  '  and  were  presented  at 

3 

the  1982  annual  meeting  of  the  Division  of  Plasma  Physics. 

They  are  briefly  reviewed  herein. 

Beam  formation  in  the  injector  has  been  studied  using 

4 

the  REEFER  code  and  the  STAGEN  code  (developed  at  Austin 
Research  Associates)  to  compute  steady  state  diode  flow 
patterns  in  r,z  geometry.  Early  examination  of  a  cantilevered 
button  field  emission  cathode  of  14  cm  radius,  2.5  cm 
thickness,  and  -400  kV  potential  inside  a  grounded  cylindrical 
housing  of  19  cm  radius  with  a  grounded  anode  mesh  located  10 
cm  downstream,  revealed  that  pre-breakdown  electric  field 
stresses  on  the  cathode  tip  would  exceed  170  kv/cm.  This  was 
found  to  lead  to  complete  breakdown  of  the  cathode  surface,  a 
current  density  profile  highly  peaked  on  edge,  some  5  kA  of 
emitted  current,  and  an  expanding  flow  pattern  in  the 
anode-cathode  gap  in  the  absence  of  a  strong  magnetic  guide 
field.  This  beam  was  judged  unacceptable  due  to  its  excessive 
size  (  >  14  cm  vs.  <  0.6  cm  desired  downstream  in  the  wiggler) , 
its  excessive  current  (5  kA  vs.  ~1  kA  desired),  and  its 
expansive  flow  in  the  diode--making  laminar  compression 


difficult  downstream.  Consequently,  the  remainder  of  our 
diode  studies  contemplated  the  use  of  either  warm  or  cold 
cathodes  with  annular  electrode  rings  in  the  anode-cathode  gap 
to  provide  a  graded  potential  axially  through  the  diode  and  to 
reduce  the  pre-breakdown  stresses  on  the  cathode  surface,  with 
emission  out  to  ~8  cm  radius  controlled  by  the  presence  of 
fibers  on  the  cathode. 

Such  controlled  emission  cathodes  were  found  to  produce 
reasonable,  high  quality  beams  with  laminar  flow  patterns. 
Diode  parameters  selected  were 

anode-cathode  gap  spacing  =  d  ®  11.0  cm 
cathode  emission  radius  *=  a  =  8.0  cm 
cathode  potential  *  -A4>  *  -400  kV. 

The  gap  spacing  of  11  cm  is  sufficient  to  minimize  the  problem 
of  diode  closure  at  velocities  of  —2.5  cm/  sec  during  the 
1-2  usee  pulse  time  of  the  beam.  The  emission  radius  of  8  cm 
is  large  enough  -to  ensure  that  the  current  density  at  the 
cathode  surface  is  less  than  4  A/cm2,  as  required  for  reliable 
operation  and  to  avoid  the  Child-Langmuir  space  charge  limit. 
However,  the  requirement  that  the  beam  should  eventually  be 
compressed  to  a  radius  of  about  0.6  cm  in  the  wiggler  solenoid 
places  a  very  high  premium  on  techniques  for  beam  convergence. 

It  is  easy  to  see  that  complete  reliance  on  magnetic 
compression  to  achieve  this  hundredfold  increase  in  the  beam 


current  density  is  hardly  feasible,  hundreds  of  gauss  would 
be  required  to  tie  the  beam  well  to  field  lines  in  the  diode 
region  and  then  compression  to  some  tens  of  kilogauss  would  be 
needed  to  reduce  the  magnetic  flux  surface  radius  to  0.6  cm  in 
the  wiggler  solenoid.  Consequently,  it  is  desirable  to 
achieve  some  beam  convergence  electrostatically  while  the 
magnetic  field  is  weak,  leaving  a  lesser  amount  of  compression 
to  be  achieved  magnetically.  Electrostatic  beam  convergence 
is  possible  in  a  graded  potential  parallel  plate  geometry. 

The  resistively  graded  electrodes  lead  to  a  radial  boundary 
condition  in  the  anode-cathode  gap  which  may  be  modeled  as 


4>  (r=b,z)  - 
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By  varying  p,  one  may  vary  the  cathode  stresses,  the  radial 
current  density  profile,  and  the  diode  flow  pattern.  The  most 
desirable  results  were  obtained  for  a  sub-Chi ld-Langmuir 
grading  of  p  =  1.52,  which  produced  relatively  weak  cathode 
stresses,  current  of  1^  e  720  A,  and  a  convergent  flow 
pattern  with  axially  peaked  current  density  as  illustrated  in 
Figures  2  and  3.  The  convergent  beam  flow  occurs  because  the 
graded  electrode  structure  tends  to  short  out  the  radial 
electric  fields  normally  associated  with  the  beam  space 
charge,  such  that  the  vz  x  B0  pinch  force  is  left  to 
compress  the  beam  radially. 
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As  shown  in  Figure  1,  the  electron  beam,  once  formed, 
must  be  transported  through  a  complex  array  of  magnetic  field 
coils  and  radially  varying  conducting  boundaries,  and  through 
two  accelerating  gaps  in  the  induction  module  before  being 
finally  compressed  and  injected  into  the  wiggler  solenoid, 
which  begins  380  cm  downstream  of  the  cathode  surface. 

In  order  to  study  this  challenging  transport  problem,  we 
have  employed  in  addition  to  theoretical  analysis,  a  one¬ 
dimensional  R(z)  beam  envelope  code  and  the  two-dimensional 
(r,z),  steady  state,  relativistic  particle  simulation  code 
STAGEN . 

The  STAGEN  code  retains  all  three  components  of 
velocity  and  of  the  electric  and  magnetic  fields,  but  is 
restricted  to  the  two  independent  variables  r,z,  appropriate 
to  situations  which  are  symmetric  in  the  azimuthal  coordinate 
6.  In  order  to  reduce  the  computing  requirements  to  a 
manageable  level,  the  axial  transport  is  divided  into  four 
"stages."  The  first  stage  extends  to  z  =  23  cm,  encompassing 
the  11  cm  length  of  the  anode-cathode  region  plus  12  cm 
beyond  the  anode  mesh.  The  flow  is  initiated  from  an 
equipotential  near  but  not  on  the  cathode  surface — typically 
4>  -  4>  ..  .  <20  kV.  Stage  2  extends  from  z  =  23  cm  to 
z  *=  111  cm,  encompassing  the  initial  focusing  coil  and  the 
conductor  wall  flaring  out  to  a  radius  of  17.5  cm.  Stage  3 
extends  from  z  =  111  cm  to  z  =  258  cm,  and  encompasses 
three  more  focusing  coils,  the  first  induction  accelerating 


gap  (of  200  kV) ,  and  the  first  induction  solenoid.  Stage  4 
extends  from  z  =  258  cm  to  z  =  406  cm,  and  encompasses  the 
second  induction  solenoid,  the  second  induction  accelerating 
gap  (also  of  200  XV),  another  focusing  coil  and  the  first 
26  cm  of  the  wiggler  solenoid. 

Under  this  contract,  STAGEN  has  been  exercised 
repeatedly  to  examine  the  beam  transport  through  the  first 
four  meters  of  the  long  pulse  induction  linac  for  a  variety  of 
diode  potential  gradings  and  for  several  magnetic  field 
configurations.  Because  the  constraints  of  operating  with  the 
existing  coil  configuration  and  achieving  an  order  of 
magnitude  beam  radius  compression  into  the  wiggler  solenoid 
make  the  transport  so  difficult,  emphasis  on  the  STAGEN  runs 
was  first  placed  on  simulating  typical  experimental 
parameters.  In  this  way,  support  for  the  experiment  could  be 
provided,  feedback  could  be  received,  and  an  appreciation  for 
the  crucial  elements  of  the  transport  physics  could  be 
obtained.  An  example  of  these  transport  simulations  is 
provided  in  Figures  4-10.  Figure  4  displays  the  magnetic 
field  profile  in  the  LPIL  for  a  typical  sert  of  coil  currents, 
including  the  magnetic  ramp  into  the  2  kG  wiggler  solenoid. 
Figure  5  shows  the  Stage  1  transport  through  the  diode. 
Theoretical  calculations  indicated  that  momentum  scatter  at 
the  anode  mesh  should  create  emittance  of  only  —25'n  mrad-cm. 
However,  experimental  indications  suggested  an  emittance  of 


200tt  mrad-cm  created  at  the  surface  of  the  hairbrush 
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cathodes.  To  mock  up  this  effect,  a  scatter  corresponding  to 

200 it  mrad-cm  emittance  was  therefore  introduced  at  the  anode 

mesh  for  these  runs.  Figure  6  shows  the  equipotential 

contours  in  the  diode  region.  Figures  7,  8  and  9  display  the 

beam  transport  through  Stages  2,  3  and  4.  In  Figure  8,  it  may 

be  seen  that  135  A  of  current  is  lost  to  the  wall  just  beyond 

the  first  induction  gap.  For  this  particular  case,  the  second 

induction  gap  was  inactive.  The  strong  beam  compression  of 

the  remaining  583  A  of  current  is  seen  in  Figure  9.  However, 

this  strong  compression  excites  large  amplitude  zero  frequency 

cyclotron  oscillations  which  are  displayed  in  the  Br~z  phase 

space  plot  of  Figure  10.  This  strong  sausaging  will  soon 

phase  mix  to  create  an  emittance  of  ~940it  mrad-cm  for  this 

example.  In  Figure  11,  the  variation  of  y  with  z  is 

z 

displayed  for  the  electrons  in  Stage  4.  At  the  throat  of  the 

wiggler,  the  median  value  of  Y  is  found  to  be  «1.6?  there 

z 

is  a  considerable  spread  since  the  outer  electrons  have  been 
excited  to  much  higher  transverse  velocities  than  have  the 
inner  electrons.  These  details  on  the  phase  space  distribu¬ 
tion  of  the  electrons  turned  out  to  be  quite  useful"*  in 
interpreting  that  certain  of  the  radiation  measurements  indeed 
signify  an  FEL  mode,  as  expected. 

The  chief  difficulty  in  this  beam  transport  is  the 
large  emittance  growth  experienced  during  the  magnetic 
compression  at  the  mouth  of  the  wiggler  solenoid.  Guidance 
for  improving  the  transport  in  this  respect  was  obtained  from 
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analysis  with  a  one-dimensional  envelope  equation  for  R(z), 


similar  to  the  equation  developed  by  Lee  and  Cooper.  The 


assumptions,  equations,  and  symbols  used  in  this  analysis  are 


summarized  in  Tables  1  and  2.  Table  3  summarizes  the  relative 


significance  of  the  various  radial  "force"  terms  in  the 


2 

envelope  equation.  The  y'  term  is  seen  to  represent  the 


focusing  kicks  that  the  beam  receives  in  the  diode  and  in  the 


induction  gaps.  During  the  early  portion  of  the  transport. 


while  B  <  200  G,  the  beam  radius  is  determined  by  competi- 
z 


tion  between  the  magnetic  v(j)Bz  focusing  force  and  the 
defocusing  space  charge  force  [terms  (1)  and  (3)  of  Table  3]. 


However,  as  the  magnetic  field  is  increased  toward  2  kG,  the 


defocusing  space  charge  force  is  outweighed  by  the  centrifugal 


and  pressure  forces  [term  (4)  of  Table  3]. 


In  Table  4,  a  Hamiltonian  analysis  of  the  cyclotron 


oscillations  excited  by  these  dominant  forces  is  summarized. 


The  beam  envelope  oscillations  may  thus  be  viewed  as 


anharmonic  oscillations  in  the  pseudopotential  P(R), 


illustrated  in  Figure  12.  This  theoretical  analysis  suggests 


that  large  cyclotron  oscillations  of  the  magnitude  observed 


should  indeed  be  created  if  the  magnetic  field  is  suddenly 


increased  greatly  at  a  point  where  the  beam  radius  is 


relatively  large.  On  the  other  hand,  it  should  theoretically 


be  possible  to  minimize  the  oscillations  by; 


$ 


Table  1 


ONE -DIMENSIONAL  ENVELOPE  EQUATION’S  FOR  TRANSPORT  OF  A 
RELATIVISTIC  ELECTRON  BEAM  THROUGH  VACUUM 

Reference:  E.  P.  Lee  end  R.  K.  Cooper.  Particle  Accelerator  7.  83(1976) 
Assumptions 


(1) 

B  f  B  (r) 

2  2 

-  neglect  diamagnetism  and  other 
radial  gradients  in 

(2) 

a 

0 

- 

-  assume  azimuthal  symmetry 

(3) 

y  t  y(r) 

- 

-  neglect  transverse  mass  spread 

(4) 

|e  1  <<  le  ! 

|  X  I  1  2  | 

- 

-  paraxial  ray  approximation 

(5) 

at 

0 

- 

-  assume  steady  state  behavior 

(6) 

J2  *  Vr):  Er 

Vr  • 

-  assumed  for  convenience 

Equations 

y(z)  • 

.  **w.ll(z) 

Jb 

2  . 

me  3  mc*/e 

6  *•*-  R(r) 

• 

•»  Y  .  * 

r  +  tr~  R  + 
8jy 

r  *„■<«>  i 

,"1 

R  Xb  _  1 

^  2 y6 2 J 

y’e’mc’/e  R  yjB3  R* 

or 

’  oo’<*) 

(yj+2)y' 

H  i  i 

ft  + 

4y*8jcj 

4y*B* 

R  - 

mc*/e  YJfiJ  fi  Rs 
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Table  2.  Definitions 


•  _  _d_ 

dz 

(Yo-l)mc2  -  electron  bees  energy  at  the  anode  potential 

6"  (l *  Y  *  6z“  normalized  axial  velocity 

R(z)  -  rms  beam  radius 

a(z)  *  Jl  R(z)  “  outer  beam  radius 

R(z)  -  (Yfi)%  R(z) 

b<z)  “  conducting  wall  radius 

•*wall(z)  •  wall  potential  (varies  in  acceleration  regions) 
nc*  ■  511  keV 


-  beam  current 


nc’/e  -  17  kA 
eBz(z) 


i  (z)  - 
o  me 


:  -  yR  «  B 


i.  rms 

n0(*) 


•  nonrelativistic  cyclotron  frequency 
-  normalized  beam  emittance  (rad-cm) ,  a  constant 


’ .  •  Y  R  B .  -  ~V  R2  •  canonical  momentum  of  beam,  a  constant 
f  9  ZC 


“  BR*  ”  normalized  radial  fluid  velocity  of  beam 


n  (z) 


B^  “  2y'c~  *  +  yR  *  normalized  azimuthal  fluid  velocity  of  beam 


[•>  c*  1% 

B2R  +  B^2  +  “  normalized  transverse  velocity 


-  normalized  thermal  velocity 
m  {Ve1***’ 3  +  Y^e^r  “  "fluid  emittance' 


Table  3 


RELATIVE  SIGNIFICANCE  OF  TERMS  IN  ENVELOPE  EQUATION 


Y2 62c!  R 


—  Magnetic  v^B2  focusing  forces  which  are 
important  everywhere  except  in  the  diode 
region  where  Bz  is  very  weak. 


(xiiiirll 

4y**  61* 


ft 


—  Radial  focusing  forces  which  are  significant 
(for  beam  current  <<  Alfven  current)  only  in 
the  diode  and  in  the  induction  gaps.  These 
forces  represent  the  replacement  of  Er  by 
Ez  fields,  and  tend  to  locally  unbalance  the 
radial  force  equilibrium. 


y2  f32mc3/e 


1 

R 


—  De focusing  space  charge  forces  which  are  in 
competition  with  magnetic  focusing  forces 
during  the  first  three  meters  of  large  radius 
beam  transport. 


c  2  +  P  2 

_ £  —  Centrifugal  and  pressure  forces  which  are  in 

competition  with  magnetic  focusing  forces 
during  the  strong  magnetic  compression  of  the 
beam  to  a  very  small  radius.  Both  c  and  P^ 
are  constant  downstream  of  the  anode,  and 
'diode  ~  200-300  mrad-cm  dominates  the 
Pg,  =  50  mrad-cm  due  to  the  weak  leakage 
magnetic  field  at  the  cathode. 
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Table  4 


ANALYSIS  OF  INCREASE  IN  ZERO  FREQUENCY  CYCLOTRON 
OSCILLATIONS  DURING  MAGNETIC  COMPRESSION 


ASSUME  Y 


V2> 

+  4YJe3c>  R  ■ 


Y3B3mc 3/e 


1 

R 


Y2  B2R3 


LET  H  *=  —  p2  +  P (R)  «  Hamiltonian 


P  (R) 


n  2<z> 

0 _ 

8y2  82c2 


R  - 


y563mc3/e 


&*iR  + 


2y2  82  R2 


R  ,  q  <-»  R,  t  *-»  z 


(1)  If  n  '(2)  “0,  H  «  constant  of  motion  “  P(R  )  •  P(R  )  s  P_ 

O  +  —  o 

R+,  R  *  turning  points  of  R(z) 

(2)  If  n  (2)  increases  adiabatically ,  then  H  «  P  (z)  increases 

O  0 

adiabatically  such  that  J  •  j  pdq  is  an  adiabatic  invariant 

j  2yBc  P0(z)  c 

When  Ifa  term  is  negligible,  -  -  - - ^ 

o 

and  one  finds  R  2  “  (z) ,  R  3  *  U  (2) 

max  o  +0 

s  n0*R* 

so  that  c/ m  YjB2R2r'  •+  »  adiabatically  invariant 

and  emittance  should  not  be  enhanced. 

(3)  If  fMz)  increases  suddenly,  the  constant  H  is  incremented 

by  fiH  *=  - — — I  hfi  2  before  R  and  R*  can  change, 
o  y  d  c  o 

Thereafter,  H  «  constant.  When  fin  2  is  large,  oscillations 

0 

are  enhanced  if  fiH!(z)  occurs  where  R2  is  large.  If 

0 

R  ■  R^  where  fil)  2  occurs,  F+  is  unchanged,  but 

R  2  «  0  J(2)  ,  so  that  e  2  _  n  <z  )2  R  Vc  2  ,  and  emittance 
max  0  f  c  *  + 


is  greatly  enhanced. 


Figure  12.  Pseudopotential  for  Beam  Envelope  Oscillations 


& 


a)  timing  a  sudden  increase  in  magnetic  field  to  occur 
when  the  beam  envelope  is  near  a  radial  minimum;  or 

b)  increasing  the  magnetic  field  strength  only 
adiabatically  in  z. 

Some  attempts  were  subsequently  made  experimentally  to 

comply  with  these  theoretical  recommendations — subject  to  con¬ 
straints  of  the  existing  coil  geometry.  An  additional  transport 

coil  was  added  to  make  the  compression  more  adiabatic  and  some 
improvement  in  transport  was  seen.  However,  the  continuing 
necessity  for  such  substantial  magnetic  compression  and  the 
constraints  of  working  with  the  existing  accelerator  configura¬ 
tion  imply  a  limitation  to  the  improvements  which  may  be 
practically  realized  in  this  way. 
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B  Beam  Emittance  Due  to  Excitation  and  Phase 

a  Nixing  of  Cyclotron  Oscillations 
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BEAM  EMITTANCE  DUE  TO  EXCITATION  AND 


PHASE  MIXING  OF  CYCLOTRON  OSCILLATIONS 


Consider  a  cold,  fluid,  relativistic  electron  beam  which 
is  magnetized  and  is  smoothly  propagating  in  a  "slowly  rotating” 
beam  equilibrium.  Suppose  that  the  electron  beam  propagates 
axially  through  a  spatial  magnetic  discontinuity.  This  will  set 
up  cyclotron  oscillations  upon  the  beam  which  will  ring  down¬ 
stream  of  the  discontinuity  until  they  phase  mix  into  an  effec¬ 
tive  beam  temperature,  or  transverse  emittance.  We  wish  to 
compute  an  expression  for  this  cyclotron  emittance  as  a  function 
of  the  magnetic  discontinuity  and  the  other  equilibrium 
parameters. 

For  the  purposes  stated  above,  we  may  employ  the  beam 
envelope  equations  to  compute  the  amplitude  of  the  cyclotron 
oscillations  prior  to  phase  mixing.  Then  an  average  of  the 
cyclotron  oscillations  over  the  axial  coordinate  may  be  performed 
analytically  to  recover  the  effect  of  phase  mixing. 

The  beam  envelope  equation  may  be  expressed  in  the  "edge 
notation"  form 


(eBz/mc) 

4y 2  B2  c2 


r  - 


2  lb  1 

y 2  B  3  me  Ve  r 


'  th 


2  +  p  2 


y2  e2 


-T  =  0 


2  v 

Y7! 


(1) 
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where 

r  (z) 

■  radius  of  the  outer 

edge  of  the  beam 

« 

_d_ 

“  dz 

I  o)  2  r2 

,  b  p 

|Ib(kA) 

V 

8mc3/e  4c2 

176 

(2) 

(eB  /me) 

P0 

-  Y‘8e  +  L 

normalized 

canonical 

momentum 

eth  -  Yrseth  *  normalized  edge  emittance 
It  may  be  shown  that 

Pg'  -  0  for  9-symmetric  beams  (3) 

*  0  for  vacuum  propagation  (4) 

and  it  has  been  assumed  that  y'  *  y"  *  0. 

The  coefficient  v  is  positive  for  electron  beams  traveling  in 
either  direction.  We  have  assumed  that  the  beam  is  traveling  in 
the  z  direction,  so  that  6  -  6  > 0.  For  electrons,  e  =  — | e | 
so  that  the  vector  potential  piece  of  P0  is  < 0  when  Bz  >  0. 

6  represents  the  normalized  9-velocity  of  the  beam 
9  as  a  whole,  and  is  characteristic  of  an  outer, 
edge  electron. 


represents  the  normalized  thermal  velocity  of 
the  beam  as  a  whole,  and  is  also  characteristic 
of  an  outer,  edge  electron. 


It  should  be  noted  that  the  above  quantities  differ  by  factors  of 


2  or  >J~2  from  the  corresponding  quantities  in  the  envelope 


equations  of  Lee  and  Cooper  (Ref.  1).  This  is  because  Lee  and 


Cooper  used  R(z)  *  rms  radius  of  the  beam,  and  correspondingly: 


r(z) 


R(z) 


30(edge) 


6e(R) 


c(z)//T 
Bq  (edge )/ yj~2 


P0(edge ) 


P0(R) 


P0 ( edge )/2 


e..  (edge) 


£th<R) 


eth(edge)/2 


The  emittance  which  appears  in  the  envelope  equations  is 


defined  in  terms  of  an  average  over  the  beam  particles  of  the 


square  of  the  transverse  particle  velocities.  For  example,  in 


the  "rms  notation"  of  Lee  and  Cooper,  one  has 


eth(R)2  s  Y2R2 


[<6/>  - 


3ZR'  2  -  B0  (R)2  =  y2R2 5 B  2 

S  rms 


- 2  d  2  x  a  2 


where 


00 

<’■> '  1 1 


dr  2"rr  J,  (r) 
b 


2  +  B  21 


i  rms 


< 


This  equation  may  be  ceexpressed  in  "edge  notation"  as  follows: 


£ 


2 

tot 


Y2r2  8 


i  edge 


where 


8  .  »  2  6 
i  edge  i  rms 

r  *  edge  radius  *  y/~2  R 


e(edge )  *  2  e ( rms ) 

dr 

Br'  *  6  r  *  6  ( r )  *  8  (edge)  *  normalized  "fluid" 

z  radial  velocity  of  the 

beam  as  a  whole, 
characteristic  of  an 
outer,  edge  electron 

6 _ ( r )  ■  normalized  "fluid"  azimuthal  velocity  of  the 

0  beam  as  a  whole,  characteristic  of  an  outer, 

edge  electron 

€th(edga) 

58th(r)  =  -  *  normalized  thermal  velocity 

Yr  of  the  beam  as  a  whole, 

characteristic  of  an  outer 
edge  electron. 


Note  that: 

edge  “  $r(edge)  +  ^g(r)  +  ^th^r^ 

The  emittance  which  appears  in  the  envelope  equation  is 
£th*  It:  character izes  the  beam  temperature  and  is  a  constant  of 
the  motion  (i.e.,  eth'  *  0)  for  the  case  of  vacuum  beam 
propagation — under  the  assumptions  used  in  deriving  the  envelope 
equations .  However,  these  assumptions  included  such  things  as 
neglecting  the  radial  dependence  of  Y  and  6  .  Therefore,  it 
is  possible  for  cyclotron  oscillations  to  be  excited  upon  the 


electron  beam,  which  are  initially  fluidlike  (i.e.,  nonthermal) — 
but  which,  because  of  the  radial  variations  in  y  and  B  , 

Z 

eventually  lead  to  phase  variations  in  the  cyclotron  motion  as 
the  beam  propagates  axially  downstream.  If  the  radial  averages 
required  in  the  emittance  definition  are  repeated  at  an  axial 
position  downstream  of  the  cyclotron  excitation,  the  corre¬ 
sponding  "phase  mixing"  will  lead  to  an  increase  in  the  beam 


I 


temperature  or  emittance,  which  is  not  predicted  strictly  from 
the  envelope  equations  alone  for  the  reason  that  these  equations 
neglect  the  radial  dependences  which  produce  the  cyclotron  phase 
variations.  However,  it  is  possible  to  compute  the  amplitude  of 
the  initial  fluid  cyclotron  oscillations  from  the  envelope 
equations,  and  this  is  sufficient  to  imply  the  post-phase-mixing 
emittance  since  the  effect  of  the  phase-mixing  averages  is  easily 
obtained  by  analytically  performing  averages  over  the  axial 
variation  in  the  cyclotron  oscillations. 

For  the  specific  problem  which  we  have  posed,  we  suppose 
that  the  electron  beam  is  initially  cold  in  the  region  upstream 
of  the  magnetic  discontinuity,  so  that  -  0.  However,  the 

magnetic  discontinuity  will  excite  fluid  cyclotron  oscillations, 
for  which  the  corresponding  fluid  emittance  (in  edge  notation)  is 

e  2  (z)  =  y2S2r2(z)r'2(z)  +y2r2(z)S2(z)  =y2r2(z)6  2  ( z) 

r  0  1  r 

(7) 
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The  eventual  downstream  thermal  emittance  which  results  from  the 


phase  mixing  of  these  fluid  oscillations  may  be  computed  by 
performing  a  z-average  to  correspond  to  the  effect  of  the  phase 
mixing,  and  then  subtracting  off  the  residual,  nonthermal  portion 
of  e^2  *  Y2r2Bl£2  which  is  due  to  the  mean  fluid  precession  of 
the  beam  in  the  downstream  equilibrium.  This  leads  to 


Aeth  “  Y282  ^r2  (z)r' 2  (z)^  +  y2<(r2  (z)  BQ2  (z)^  -  Y2  <^r  (z)  BQ  (z)^> 

z  z  z 

(8) 


1  2 

where  y2  (z )  B0(z  is  the  component  of  £f2  *  Y2r2Bif2  due 

to  the  mean  fluid  precession  of  the  beam  in  the  downstream 
equilibrium.  Here  represents  a  z-average  over  one  cycle 

)  of  the  cyclotron  oscillations. 

j  By  employing  the  canonical  momentum  (constant) 

i 

I  Equation  (2),  one  may  eliminate  Bg(z)  in  favor  of  r(z),  so 

that  the  increase  in  thermal  emittance  may  finally  be  expressed 
in  terms  of  the  beam  radius  function  r(z)  alone. 


I 


n 


Y  r (z)  Bq (z) 


0  o 

P  +  —  r2  Cz) 
0  2c 


(9) 


e  B 


where 


a 


o  me 

It  follows  that 


>  0  . 


aSh  =  + 

z 


a 


0 

4c2 


[<r>(z)>  -  <rMz)>z] 


(10) 
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It  remains  only  to  compute  the  envelope  oscillations 


which  are  excited  by  a  magnetic  discontinuity,  and  then  to 
evaluate  Ae^2  which  results  from  the  phase  mixing  of  these 
oscillations. 

Roberson  (Ref.  2)  linearized  the  envelope  Equation  (1) 
and  solved  for  the  linearized  envelope  response  to  the  magnetic 
perturbation  shown  in  Figure  1,  which  is 


Figure  1.  Magnetic  discontinuity  which  excites 

cyclotron  oscillations  upon  an  electron  beam 


ft  (z) 

o 


I e |  B  (z) 

z 

m  c 


ft  ,  z  <  0 

o 

ft0  +  6ft,  0  <  z  <  L 
ft  ,  L  <  z 


(1 


|  e  |  6  B 

z 

5ft  =  - 


where 


and  6ft  can  be  of  either  sign.  The  solution  to  Equation  (1)  is 
found  to  be 


;  z  <  0 


r  (z)  =  <  r  1 


[*- 


ft  6ft 
o 

,  o  9".  j  o  sin 
k2Y2B2c2 


<  z  <  L 


ft  6ft 
0 


k2Y262c2 


(¥)]  -  0 

kL  .  /.  kL\l 

Ln  t  sinlkz  ’  T/J  ; 


where 


ft 


y282c2 


4v 

Y  382rQ2 


In  the  downstream  region  L<z, 
oscillations  have  the  form 

r(z)  =  rQ  +  6r  sin  |kz 


the  linearized  cyclotron 


where 

6r 


r  ft  5ft 
0  0 


k2Y282c2 


sin 


kL 


If  this  form  is  inserted  into  the  expression  (10)  for  Ae 
the  result  is  evaluated  to  lowest  order  in  5r2/rQ2  ,  one 


(12) 

i  <  z 

(13) 

(14) 

(15) 

th 

obtains 


2 


<6r2 (z)> 
z 

(16) 


As  an  aside,  consider  for  a  moment  the  beam  equilibrium.  If  the 
linearized  form  r  *  rQ  +  6r  sin(kz  -  kL/2)  is  inserted  into  the 
envelope  equation  and  the  equation  is  then  averaged  in  z,  there 
results 


where  the  equation  is  carried  to  9\6r2/zQ2  )  to  suggest  the 
influence  of  <5r2  upon  the  equilibrium.  The  assertion  that  the 
RHS  above  must  be  positive  gives  one  of  the  equilibrium  criteria. 
In  the  limit  that  5r  »  0  and  efch  *  0,  this  reduces  to  one  of 
the  usual  criteria  for  cold,  fluid  equilibria: 


The  second  equilibrium  criteria  comes  from  requiring  that  the 
mean  precessional  velocity  should  obey  ]30  <  yj  1  -  -  1/y 

as  shown  above.  It  maybe  seen  that  the  presence  of  waves  (i.e., 
<Sr 2  f  0)  causes  the  mean  precessional  velocity  to  increase,  and 
hence  makes  this  second  equilibrium  criteria  more  difficult  to 


satisfy  as  well.  In  the  limit  of  6r 


£ th •  this  criterion 


reduces  to  the  second  criterion  for  cold,  fluid  equilibria: 


,  2v  c 

1  *  T  — 

1  o  o 


rb(kA) 

5  B  (kG)  r  (cm) yg 
z  o 


and  for  large  values  of  fir /c,  this  criterion  is  the  more 


demanding. 


If  the  above  expression  for  is  employed,  one  may 


show  that 


<vr(z)  !!„(»))* 

2 r 2  ( z )  8q2(z)^>z 


.2  ,  2vc“ 


e.  2c2 


Sr  +  K  !  r  s 

0  0  0 


o  e  -  c 

1  ?  -  2  .  2vc2  .  th 

—  r  or  +  — — 7  +  — — 2 - 5 

2  o  y  ft  2  ft  2  r  2 

0  0  0 


where  the  6r2/r  2  corrections  to  v  and  e  2  have  been 

o  tn 

omitted. 

When  6r2  -  0,  <yr  80>2  —  <Y2r2B02>z  so  that 

2  2  2  /  v  ^ 

Y  r  8g  is  dominated  by  the  fluid  piece  <Qyr60^>  ,  and 


Ac..  "*  0*  However,  if  the  beam  is  excited  to  the  degree  that 


g 


then 


<^2V> 


(24) 


i 


(25) 


so  that  the  fluid  precessional  piece  of  y2r2£3  2  is  negligible 
in  comparison  with  the  cyclotron  wave  piece. 

When  the  magnetic  guide  field  is  strong  enough  that  the 
equilibrium  criteria  are  very  well  obeyed,  the  increased  thermal 
emittance  due  to  the  phase  mixing  of  cyclotron  oscillations  is 
given  by 


A  e 


th 


2fl  2  r  2 
o  o 

- — T5 — 


<5r2(z)>  = 


5r 2 


(26) 


If  this  "cyclotron  wave  emittance"  is  the  dominant  source 
of  parallel  velocity  spread,  then 


y26  A8 

z  z 


aeth2  no2<5r2(z)>z 

2  r  2  c2 
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A  BRIEF  NOTE  ON  BEAM  QUALITY  IN  HIGH 
CURRENT,  FIELD  IMMERSED  DIODES* 

M.  L.  Sloan  and  J.  R.  Thompson 

I .  INTRODUCTION 

Analytical  calculations  and  numerical  particle 
simulations  have  indicated  that  very  significant  improvements 
in  the  quality  of  high  current  relativistic  electron  beams 
are  possible  through  the  use  of  appropriately  shaped  diode 
structures  which  are  immersed  in  very  strong,  longitudinal 
magnetic  guide  fields.  Improvements  in  the  beam  emittance 
of  one  to  two  orders  of  magnitude  over  the  Lawson-Penner 
emittance 

en  =  y6R0  =  320  Ifa(kA)  mrad-cm  (1) 

have  been  observed  in  some  cases. 

Increasing  the  strength  of  the  axial  magnetic  field  can 
in  itself  reduce  the  transverse  velocity  acquired  by  beam 

*This  work  was  supported  by  the  Office  of  Naval  Research 
under  Contract  No.  N00014-81-C-0704. 


electrons  in  traversing  a  given  diode  structure.  As  shown  in 
Section  II,  theoretical  estimates  are  that  magnetic  fields 
sufficient  for 


1  << 


eB  Az 
z 


yB  mc‘ 
z 


are  required  for  significant  beam  quality  improvement,  where 


Az  is  the  gradient  length  of  transverse  diode  forces. 


Depending  upon  the  type  of  diode  structure,  the  transverse 


electron  velocity  perturbations  acquired  may  scale  as 


(ftAz/B  c)-i  or  as  exp(-nAz/8  c) . 
z  z 


For  a  given  magnetic  guide  field  strength,  sufficient 


that  the  electrons  are  well  tied  to  field  lines,  beam  quality 


may  be  strongly  enhanced  or  degraded  by  changes  in  the  shape 


of  the  anode  and  cathode  structure.  The  effects  of  such 


changes  may  be  explored  via  numerical  particle  simulations. 


Austin  Research  Associates  has  developed  a  code  DIODESIGN 


which  is  very  useful  for  such  high  magnetic  field  diode 


studies,  since  it  allows  a  direct  calculation  of  the  shape  of 


the  cathode  structure  required  to  produce  a  desired  current 


density  profile.  It  is  generally  found  that  the  best  beam 


quality  results  from  gently-shaped  anode  and  cathode 


LM .  L.  Sloan  and  H.  A.  Davis,  Phys.  Fluids  25^  2337  (1982) 
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structures  which  create  electric  fields  whose  radial 

components  have  axial  gradient  lengths  long  compared  to  the 

inverse  cyclotron  wave  number  Bzc/Q. 

To  illustrate  these  ideas,  we  have  performed 

several  calculations  for  a  diode  structure  similar  to 

the  VEBA  apertured  diode  which  has  been  employed  in  free 

2 

electron  laser  experiments  at  NRL.  Figure  3,  taken  from 
Reference  2,  illustrates  the  VEBA  geometry.  It  is  reported 
that  with  this  geometry,  with  the  guide  field  strength  set  at 
15  kG,  6a  ~  0.034  for  the  transmitted  beam  electrons. 

Figure  1  displays  the  geometry  adopted  for  our  computer 
simulation  studies.  In  Table  1,  parameters  are  listed  for 
three  simulations  which  have  been  selected  to  illustrate  the 
effects  of  magnetic  field  strength  and  of  diode  geometry  upon 
the  transverse  velocity  acquired  by  the  beam  electrons. 

It  may  be  observed  that  for  a  fixed  smooth  anode  taper 
(i.e.,  e-folding  length  =s  0.51  cm),  increasing  the  magnetic 
field  from  7.44  kG  to  15.0  kG  results  in  the  average 
decreasing  from  0.0221  to  0.0024 — a  factor  of  nine!  In  this 
diode  geometry,  it  is  believed  that  the  primary  source  of 
transverse  velocity  perturbations  is  the  axially  varying 
radial  electric  field  near  the  entrance  to  the  drift  tube. 

As  shown  in  Section  II,  such  downstream  perturbations  should 

2Robert  H.  Jackson,  et.  al.,  IEEE  J.  Quantum  Electron.  QE-19 , 
346  (1983)  . 
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TABLE  1 


PARAMETERS  FOR  BEAM  QUALITY  STUDIES 
IN  A  MAGNETIZED,  FOILLESS  DIODE 


Low  Bz, 

Smooth  Anode  Taper 


High  Bz , 

Smooth  Anode  Taper 


High  Bz, 

Sharp  Anode  Taper 


1.5  MV 

1.5  MV 

1.5  MV 

magnetron 

1.81  kA 

0.47  kA 

0.47  kA 

^sheath 

6.3  kA 

8.0  kA 

8.0  kA 

B 

1.2  kA 

1.2  kA 

1.2  kA 

5.20  cm 

5.20  cm 

5.20  cm 

o 

0.85  cm 

0.85  cm 

0.85  cm 

f 

0.306  cm 

0.306  cm 

0.306  cm 

tot 

15.20  cm 

15.20  cm 

15.20  cm 

drop 

6.80  cm 

6.80  cm 

6.80  cm 

c 

6.60  cm 

6.60  cm 

5.64  cm 

0.51  cm 

0.51  cm 

0.034  cm 

z 

7  44  kG 

15.0  kG 

15.0  kG 

0.0221 

0.0024 

0.0605 

<b.>  af 

27  mrad-cm 

3  mrad-cm 

73  mrad-cm 

0.87  cm 


0.43  cm 


0.43  cm 


fall  off  roughly  exponentially  as  the  magnetic  field  strength 
is  raised. 

It  may  also  be  observed  that  for  a  fixed,  strong 
magnetic  field  of  15  kG,  a  dramatic  sharpening  in  the  anode 
taper  length  (from  0.51  cm  to  0.034  cm)  results  in  the 
average  8^  increasing  from  0.0024  to  0.0605 — a  factor  of 
twenty-five!  What  has  happened,  of  course,  is  that  a  sharp 
corner  has  been  created  near  the  mouth  of  the  drift  tube, 
which  creates  radial  forces  with  a  short  axial  gradient 
length.  The  sharpness  of  the  corner  in  this  third  simulation 
is  actually  a  bit  worse  than  the  actual  corner  shown  in 
Figure  3  for  the  VEBA  diode,  which  accounts  for  the  simula¬ 
tion  6a  of  0.0605  being  almost  twice  the  value  reported  in 
the  actual  VEBA  experiment.  Nevertheless,  it  seems  clear 
that  the  experimental  beam  quality  could  undoubtedly  be 
improved  significantly  if  the  anode  taper  were  made  more 
gentle. 

Our  studies  suggest  that  by  employing  magnetic  fields 
strong  enough  that  beam  electrons  are  tied  to  the  field  lines 
and  that  Bzc/ft  is  short  compared  to  the  axial  gradient 
length  of  radial  forces,  and  by  designing  smooth  anode- 
cathode  structures,  extremely  high  quality  electron  beams  may 
be  created.  Although  it  has  not  been  considered  in  this 
note,  it  should  also  be  recognized  that  before  such  high 
quality  beams  can  be  exploited,  it  is  necessary  to  design  a 


beam  transport  system  capable  of  transporting  the  beam  to  the 
point  of  use  without  degrading  the  previously-acquired 
quality. 


II.  DEPENDENCE  OF  BEAM  TEMPERATURE  ON 


MAGNETIC  GUIDE  FIELD  STRENGTH 

i 

l 

t 

! 

We  assume  an  axisymmetric  diode  operated  in  a  strong  j 

i 

longitudinal  magnetic  field  B  -  BQz  where  the  beam  is  j 

I 

sufficiently  cold  (or  conversely  BQ  sufficiently  large) 

that  the  Larmor  radius  («p  /B  )  is  small  compared  to  the 

J.  o 

dimensions  or  gradient  lengths  of  the  system.  We  also  assume  j 

that  to  lowest  order  the  slow  beam  rotation  is  negligible  so 

that 

£  Pzz 


dz  _ 
dx 


pz  = 


Y6, 


d  _  d 
dx  =  Y  dt 


(3) 


where  x  is  the  proper  time  measured  along  the  beam  orbit. 
For  such  a  strong,  straight  magnetic  field,  the  beam  electron 
radial  displacement  can  be  neglected  to  lowest  order,  and  the 
beam  velocity  perturbations  may  be  computed  by  integrating 
the  force  equations  along  straight  line  orbits.  The  perpen¬ 
dicular  force  equations  are 


(4) 


-  p.ft  +  F 
$  o  r 


p  ft  +  F 

r  o  <p 


where  ft  =  eB„/mc,  e  =  lei  ,  and  F  ,  F,  are  the  normalized 

o  o  ii  r  0 

forces  on  the  electrons  due  to  the  applied  and  self  fields. 

Under  the  assumption  that  the  gyroradius  is  small  and 
the  electron  moves  in  z  according  to  Equation  (3),  we 
i ntroduce 


P"  5  Pr  1  ^ 


to  obtain  a  compact  form  for  Equation  (4) 


d  ±  ± 

—  p  ;  i!!  p  =  F  ±  IF 
dx  o  r  <J) 


which  has  an  immediate  integral 


l 

,!  e  IiI!oT  -  j dx'  (Pr  ± 


X  ift  T  ’ 
lF  J  e  o 
4> 


where  the  integration  is  taken  along  the  zero  order  path 


r  ■  constant 

dz  (7) 

dx  * 

and  the  origin  t  *  o  is  taken  at  the  cathode  surface. 

The  initial  condition  p±|t=q  =0  is  assumed. 

We  assume  that  far  downstream  of  the  diode  region 
the  beam  is  again  in  force  equilibrium,  for  which 
F  ±  iF .  -*■  constant,  beyond  some  t  ,  and  hence  for 
T  >  Tmax  the  upper  limit  of  the  integral  in  Equation  (6) 
may  be  extended  to  infinity.  Thus  the  solution  for  the 
perturbed  beam  velocities  far  downstream  of  the  diode  may  be 
written  as 


00 

±  ±ifl  T  f 

)  =  e  o  J 


.  , .  + in  t ' 

dt '  (F  ±  lFJ  e  o 
r  <(> 


Defining  p,2  =  p  2  +  pj,  we 

•L  IT  U) 


obtain 


(8) 


PA2  =  P+  P-  = 


00 

JdT  '  (Fr 


-it!  T' 

+ iF.)  e  o 


(9) 
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Thus  the  transverse  beam  temperature  is  related  to  the 
Fourier  transform  of  the  forces  acting  upon  the  electrons  in 
the  diode  region,  evaluated  at  the  gyrof requency . 


The  transverse  forces  are  given  by 


F  =  —  <-E  +6  B.) 
r  m  r  z  <f> 


t-t  «H|.| 


(-E  +  6  B  ) 
r  m  r  z  <J> 


_  PrP<|> 

F ,  =  -  - —  as 

<J>  r 


(10) 


where  the  centrifugal  and  centripetal  forces  may  be  neglected 
in  the  strong  magnetic  field.  Hence 


y\  -  - 


a 
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dt  F  (t) 
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-  ift  x 
e  o 


e 

me 


(11) 


It  is  apparent  from  Equation  (11)  that  large  transverse 
velocity  perturbations  may  result  when  the  fields  Er  -BzB^ 
have  significant  components  on  the  same  scale  length  as 
y6zc/^0.  Likewise,  increasing  Bq  should  significantly 
reduce  by  reducing  the  latter  gyrolength  below  the 

transverse  force  gradient  lengths. 


Qualitatively,  the  forces  Fr  tend  to  scale  as  the 
beam  current,  although  they  also  reflect  the  applied  diode 
fields.  These  forces  will  tend  to  be  relatively  larger  for 
the  outer  beam  electrons  than  for  axial  electrons.  If 
unused,  to-be-apertured  beam  electrons  are  present,  these 
will  tend  to  shield  the  interior  electrons  somewhat  from  the 
effects  of  applied  transverse  diode  fields. 

To  proceed  further  with  this  calculation  requires 

detailed  knowledge  of  the  diode  field  structure.  However,  we 

may  examine  two  model  calculations  which  serve  to  illustrate 

the  scalings.  The  integrand  in  Equation  (11)  may  be  seen  to 

have  variations  near  the  cathode  due  to  variations  in  8z  as 

the  electrons  are  first  accelerated,  with  the  integrand 

declining  as  8  increases.  These  cathode  effects  may  be 
z 

roughly  modeled  as 

Fr(T)  *  Fq  e  “t/Ato  (12) 

which  leads  immediately  to  the  result 

F  Ax 
o  o 

^iO  yc  (1  +  2  Ax2 )  ^ 

0  0 
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^  v 


r  r^  Wv,'*1  *  V  v 


where 


F  =  —  (-E  +6  BJ 
o  in  r  z  i> 


At  «  — - — 
o  yB  c 
z 


and  Az  is  the  acceleration  scale  length.  Since  E  is 
o  r 

normally  small  near  the  cathode,  these  8io  perturbations 
should  be  small  and  should  fall  off  as  B0“^  as  BQ  is  made 
very  large. 

The  integrand  in  Equation  (.11)  may  also  acquire 
significant  variations  downstream  of  the  cathode  near  any 
apertures  or  electrical  constrictions  near  the  beam  path. 

The  amplitude  of  these  variations  may  be  much  larger.  Such 
downstream  effects  may  be  modeled  as 


r  ^  T  1  +  (t-t^  J2/^2 


When  the  width  Atj  is  much  less  than  the  location  x  ,  the 
lower  limit  on  the  integral  in  Equation  (111  may  be  extended 
from  0  to  and  the  immediate  result  is 


■.  "v  -r-  V*  *r  r 


r  k1-,  n'. 


r  i  *  9  • t  i 
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ir  f  Atj  e  o  i 

yc 


where 


=  —  (-  E  +  8  B  ) 
m  r  z  <t> 


T  =  T. 


Az 

At  -  -■ 

i  y8  c 


(15) 


and  AZj  is  the  effective  axial  scale  length  of  the  radial 
forces  near  the  downstream  constriction  at  t * t  .  Although 
Fj  may  be  rather  large,  in  this  case  the  8X ^  perturbations 
may  be  made  exponentially  small  when  YBzc/ftQ  is  reduced 


below  Az 
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